NUSC  Technical  Report  8077 
15]  August  1987 


LIBRARY 

Ksearch  reports  nmm 

n?Jfhl  POSTGRADUATE  SCHOd 
MOHTEREY.  CA!  IFORNIA  939-20 


Space  Time  Processing, 
Environmental-Acoustic  Effects 


William  M,  Carey 

NUSC,  New  London,  CT  06320 

William  B.JMosele^ 

NORDA,  NSTL,  MS  39529 


Naval  Underwater  Systems  Center 

Newport,  Rhode  Island  /  New  London,  Connecticut 


Approved  for  public  release;  distribution  is  unlimited. 


PREFACE 


This  document  is  an  expanded  version  of  an  invited  paper  presented  at  the 
12th  International  Congress  on  Acoustics  held  in  Halifax,  Nova  Scotia,  in  July, 
1986.  The  sponsors  of  this  work  are  numerous  and  the  work  itself  spanned  the 
activities  of  several  institutions,  most  notably  the  Naval  Research  Laboratory 
and  the  Naval  Oceanographic  Research  and  Development  Activity.  The  results 
presented  in  this  document  were  necessarily  drawn  from  the  NRL  research  group 
which  was  under  the  supervision  of  Dr.  Moseley  and  numerous  contractors  under 
the  supervision  of  Dr.  Carey.  This  report,  prepared  by  W.  Carey,  is  issued  to 
document  and  to  make  available  the  results  of  several  years  of  research. 
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INTRODUCTION 


This  paper  treats  the  subject  of  space-time  processing  and  the  application 

of  ocean-acoustic  arrays.  Space-time  processing  has  been  extensively  studied 

from  a  radar  viewpoint  since  the  forties  and  many  comprehensive  reviews  and 

texts  can  be  found.  The  Winder  and  Loda  (1962)^  summary  of  acoustic-space- time 

processing  together  with  the  fundamentals  of  statistical  communication  the- 
2-5 

ory  provides  a  basis  for  the  processing  of  ocean  acoustic  waves.  The  prob¬ 
lem  of  the  response  of  the  acoustic  antenna  to  a  signal  with  random  parameters 

6  7  8 

was  first  treated  by  Bourret  followed  by  Berman  and  Bordelon  .  Additional 

9  10 

work  was  performed  by  Brown  and  Lord  and  Murphy  (also  see  ref  11).  Shif- 
12 

rin  applied  a  similar  approach  to  determine  the  statistical  characteristics 

of  radar  antennas  and  his  text  treats  in  detail  the  response  of  antenna  systems 

to  waves  with  random  amplitude  and  phase  components.  His  results  will  be  used 

13-17 

as  they  apply  to  ocean  acoustic  arrays.  Additional  reviews  and  texts 

stress  the  digital  signal  processing  aspects  of  the  problem.  These  include  the 
use  of  digital  Fast  Fourier  Transforms  (FFT) ,  shading  characteristics^5^  and 
high  resolution  techniques.^  ^  These  reviews  for  the  most  part  address  spa¬ 
tial  and  temporal  processing  separately  but  have  direct  application  to  sonar 
processing  in  space-time. 


With  such  a  formidable  amount  of  previous  papers,  reviews,  and  texts,  one 
may  question  the  necessity  for  an  additional  paper  on  the  subject.  However, 
in  underwater  acoustics,  the  changes  caused  by  multipath  propagation  and  rela¬ 
tive  source-receiver  motion  through  the  severe  acoustic  interference  field 
merit  discussion.  This  paper  utilizes  the  information  available  in  the  pre¬ 
vious  reviews  as  applied  to  the  characteristics  of  acoustic  antennas  in 

the  mid- frequency  range.  The  paper  addresses  aspects  of  the  problem  unique  to 
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the  acoustic  antenna  which  result  from  the  ocean  environment  and  the  acoustic 
field.  It  also  treats  the  problems  of  measurement  of  the  spatial  coherence 
and  the  characteristic  of  linear  arrays.  Finally  the  paper  discusses  space- 
time  processing  when  the  temporal  and  spatial  variables  are  weakly  coupled. 
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THEORETICAL  BACKGROUND 

This  paper  addresses  the  space-time  processing  of  acoustic  signals  which 
satisfy  the  wave  equation  that  in  the  absence  of  surfaces  which  scatter  or 
diffract  can  be  represented  by: 

□  2  't'Cx.t)  =  (x  ,t  ) 

—  s  — o  o 

’FCx.t)  =  ffff  dt  dV  G(x,x  ,t,t  )  4>  (x  ,t  )  (1) 

—  J  J  J  J  o  o  — — o  o  s  o 

□2  G(x,x  ,t,t  )  =  -4tt6(x-x  )  6(t-t  ) 

- o  o  - o  o 

These  solutions  (¥,G)  are  known  to  be  continuous  functions  with  continuous 
first  derivatives.  The  solutions  are  also  bounded  by  the  imposition  of  a  radi¬ 
ation  condition  at  infinity.  The  solutions  and  the  source  function,  0  (x  ,  t  ), 

s  o  o 

usually  have  the  property  that  the  temporal  and  spatial  characteristics  are 
separable.  Furthermore,  the  sources  of  sound  of  interest  here  are  harmonic  and 
have  been  excited  for  a  time  period  sufficiently  long  that  a  steady  state  has 
been  achieved.  Since  the  acoustic  pressure  can  be  written  as 

P(_x,t)  =  -p  3Y(x,t)/3t  (2) 

and  since  superposition  applies,  the  Dirichlet  conditions  are  satisfied  and 
Fourier  analysis  can  be  used  to  describe  the  space-time  properties. 

THE  FOURIER  TRANSFORM  IN  SPACE  AND  TIME 

1 8 

The  pressure  can  be  written  as: 

P(_x,t)  =  {///  Fp(w,k.)  exp(  iwt-ilc’x)  dwd3k  (3) 

and  the  frequency  -  wavenumber  transform  as: 

Fp(oj,jO  =  jjjj  p(x_,t)  exp(-ia)t+ilc*3c)  dtd3x  (4) 
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This  frequency-wavenumber  transform  contains  the  spatial  and  frequency 
characteristics  of  the  acoustic  pressure  field.  This  transform  is  related  to 
Lhe  space-time  autocorrelation  and  correlation  function  of  the  acoustic  pres¬ 
sure.  The  space- time  correlation  function  is: 


Rp(tljt2,Xi.x2)  =  ////  P(xlSt)  P  (x2,C2)  dt2d-2 


(5) 


In  the  cases  of  a  harmonic  field  which  is  steady  or  for  a  random  field 
which  is  spatially  homogeneous  and  temporally  stationary,  one  can  infer  the 
space-time  correlation  function  is  only  dependent  on  the  relative  separations 
in  space  and  time. 


R  (t1,toJx  ,x_)  =  R  (t  ,t  4r,x  ,x-+r)  *  R  (x,r)  (6) 

P  I  Z  — i  — z  p  1  1  — o  — I  —  p  — 

Furthermore,  if  an  arbitrary  functions  is  amenable  to  Fourier  analysis, 
then  according  to  the  Wiener-Khintchine  theorem  the  correlation  function  and 
the  power  spectra  are  Fourier  transforms. ^  ^  The  frequency-wavenumber  spectra 
is 

Sp(ai,k_)  =  JJJJ  Rp(x,_r)  exp[-i(u>T-k#r)  ]  dxd3_r  (7) 

and  the  correlation  function 


Rp (T  >£)  =  J//J  Sp(u»]i)  exp[i(wt-ik-r)  ]  dwd3k  (8) 

These  relationships  define  the  acoustic-space- time  field  for  the  class  of 
harmonic  and  random  functions  which  are  spatially  homogeneous  and  temporally 
stationary.  These  relationships  show  that  when  a  narrowband  characteristic 
exists  in  the  frequency-wavenumber  domain,  i.e.,  a  peaked  distribution  in  w  and 
k,  then  the  correlation  function  is  broad  in  time  t  and  with  respect  to  the 
spatial  coordinate  r.  In  this  case  the  field  is  termed  spatially  and  temporal- 
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ly  coherent.  On  the  other  hand,  if  the  temporal  correlation  function  is 
strongly  peaked  in  _r  and  t,  then  the  frequency-wavenumber  spectra  is  broad  and 
the  field  is  termed  incoherent. 

THE  STATISTICAL  REPRESENTATION 


The  usual  method  of  determining  the  correlation  function  has  been  to  use 
temporal  averaging  of  the  pressure  field.  An  alternative  method  would  be  sta¬ 
tistical  averaging  or  the  computation  of  the  expectation  value.  The  temporal 
representation  is: 

R  .(x,r)  =  lim  J-t  P+(t+r,jr)  dt  =  <  P(t,£)  P+(t+T,£)  >  (9) 

P  T-H» 

and  the  statistical  determination  is: 


Rps(t,r)  =  E  (P1(t,r)  P2(t+r,r)}  =  f j _m  P1P2P(P]L,P2)  d?1 


dP, 


(10) 


where  P (P^ jP^)  is  the  joint  probability  density  function  for  the  ensemble  of 

values  P^,  P^  and  the  process  is  a  stationary  one.  The  statistical  computation 

R  (t ,r)  is  only  equal  to  the  temporal  computation  R  (x,r)  when  the  process  is 
p  s  p  t  — 

deterministic  or  ergodic.  This  is  an  important  consideration  in  the  comparison 
of  space-time  processes  which  are  statistically  computed  and  which  are  tempo¬ 
rally  measured.  The  statistical  averaging  corresponds  to  all  the  different 
states  encompassed  by  the  ensemble;  however,  the  temporal  averaging  only  corre¬ 
sponds  to  those  states  actually  observed.  Similar  reasoning  applies  to  the 
averaging  over  the  spatial  variable,  jr.  When  the  field  is  homogeneous  and  sta¬ 
tionary,  then  (in  large  average  limits)  spatial  and  temporal  average  values 
approach  the  statistically  determined  expectation  values. 
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THE  SPACE-TIME  TRANSFORM  SUMMARY 

The  relationships  between  the  space-time  transforms  has  been  discussed 
previously  and  can  be  summarized  as  follows: 

THE  SPACE  TIME  TRANSFORMS 

F(P(r , t) ) 


P(r,t) 


V  V 


- - ►  F  (u),k) 

_  P  — 

f_1(f  U,k)) 

P  — 

?2^r2 ’ t+T  >> 

F  (oj,k)  * 
P  “ 

*  F+(u>,k) 
P  ~ 

F(R  (t , r) ) 

P  ~ 

r 

.  ►" 

S 

(w  5  k) 

F_1(S  (u ,k) ) 

P  “ 

P 

Rp(x ,r) 


where  F  and  F  ^  refer  to  the  Fourier  transform  and  its  inverse.  These  trans¬ 
formations  are  useful  in  the  determination  of  the  properties  of  the  particular 
field  P(_r,t)  provided  one  can  measure  S^(o),]c). 

To  describe  the  field  P(r,t)  by  the  measurement  of  S  (w,k)  one  needs  to 

—  p  — 

account  for  the  response  of  the  measurement  system.  If  this  response  is  desig¬ 
nated  by  f(_r,t)  which  describes  the  spatial  and  temporal  response  of  the  system 

and  F  (o),k)  represents  the  measured  transform 
mp  — 

F  (to  j  k )  =  F{  f  ( r ,  t )  P(  r ,  t )  }  =  F{f(r,t)}  *  F{P(r,t)}  =  f,  *  F  (11) 

mp  —  —  —  —  lp  p 


where  *  refers  to  convolution. 


Thus  the  frequency-wavenumber  spectra  is  the  convolution  of  the  transform 
of  the  system  response  and  the  transform  of  the  field.  This  is  a  very  useful 
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concept  as  it  allows  one  to  consider  the  system  response  and  its  Fourier  trans¬ 
form  independent  of  the  acoustic  field  variables.  Furthermore,  since 


S  (o),k)  =  F  (a),k)  F  (w,k) 
mp  —  mp  —  mp 


(12) 


represents  the  output  of  the  system,  one  concludes  that  broadness  in  the  system 

response  function  f^^(a),k)  and  broadness  in  acoustic  wavenumber  spectra  both 

determine  the  broadness  of  S^Ca^lc).  This  concept  implies  that  in  the  spatial 

domain,  the  broadness  of  f  (u),k_)  with  respect  to  Ic  is  determined  by  the  shad- 

xp 

ing  characteristic  and  length  of  the  system.  Matched  processing  would  require 


that  the  spread  of  k  for  both  f  (ai,k)  and  F  (u),k)  would  be  equal.  Recall  that 

—  *-P  P 

the  angular  spread  in  F^(co,^k)  is  caused  by  the  environment. 


SYSTEM  RESPONSE 


The  description  of  the  system  response  can  easily  be  developed  by  the  di¬ 
rect  analogy  between  sampling  in  time  and  sampling  in  space.  Since  the  devel¬ 
opments  elsewhere  concentrate  on  the  time  domain,  this  paper  will  proceed  with 
a  spatial  description  of  the  system  response.  The  particular  case  of  a  one¬ 
dimensional  linear  sample,  a  line  array,  is  described  by  the  illumination  func¬ 
tion,  shading,  or  window  function.  Several  representative  window  functions  are 
shown  in  Figure  1  based  on  the  integral  form  of  the  Fourier  transform.  The 
rectangular  window  transform  follows  directly  from  equation  4.  Neglecting  the 
time-frequency  transform,  one  has 


=  i*-a/2  exp(ikx)  dx  =  a  sin  (ak/2)/ak/2 


(13) 


the  result  is  the  same  whether  the  sampling  is  in  space  or  in  time.  For  the 
case  of  discrete  elements  within  the  array  we  have 
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Figure  1.  Window  functions  determined  from  the  Fourier  integral  for  finite  apertures. 
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N  sin(Nk  dsinG) 

Vk)  =  f  l  6(x  •  nd)exp(ikx)dx  =  Nsin(k°dsin0)  (U) 

n=l  o 

This  result  agrees  closely  with  the  continuous  line  array  (equation  13) 

when  the  spacing  is  half  wavelength  (d  =  X/2)  as  shown  in  Figure  2.^  The  win- 

14 

dow  factors  listed  in  Table  1  are  taken  from  Harris's  concise  review  of  dis¬ 
crete  Fourier  transforms  and  the  use  of  windows.  For  the  purpose  of  this  paper 
we  will  continue  with  the  integral  (continuous)  description  of  the  transform, 
recognizing  that  small  differences  do  exist. 


t 

LU 

cc 

=> 

GO 


Figure  2.  A  comparison  of  the  directivity  patterns  of  (a)  a  discrete  linear 

array  having  elements  spaced  at  half  wavelengths  and  (b)  a  continu¬ 
ous,  uniform  line  transducer. 
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Table  1. 


Representative  window  function  in  properties  after  Harris 


( 1978) 14 


WINDOW 

HIGHEST 

SIDEL0BE 

LEVEL 

(dB) 

SIDEL0BE 
FALL0FF 
(dB  /  OCT) 

COHERENT 

GAIN 

EQUIV. 

NOISE. 

BW 

(BINS) 

3.0-dB 

BW 

(BINS) 

Dirichlet.  rectangle 

-13 

-6 

1.00 

1.00 

0.89 

Fejer  triangle 

-27 

-12 

0.50 

1.33 

1.28 

cos  a  [X] 

a  =  1.0 

-23 

-12 

0.64 

1.23 

1  20 

Hann 

a  =  2.0 

-32 

-18 

0.50 

1.50 

1.44 

a  =  3.0 

-39 

-24 

0.42 

1.73 

1.66 

a  =  40 

-47 

-30 

0.38 

1.94 

1.86 

Hamming 

-43 

-6 

0.54 

1.36 

1.30 

Tukey 

-14 

-18 

0.88 

1.10 

1.01 

(cosine 

-15 

-18 

0.75 

1.22 

1.15 

tapered) 

-19 

-18 

0.63 

1.36 

1.31 

exp(-aX) 

a  =  2.0 

-19 

-6 

0.44 

1.30 

1.21 

Poisson 

Q  =  3  0 

-24 

-6 

0.32 

1.65 

1.45 

a  =  4.0 

-31 

-6 

0.25 

2.08 

1.75 

Hann- 

a  =  0.5 

-35 

-18 

0.43 

1.61 

1.54 

Poisson 

a  =  1.0 

-39 

-18 

0.38 

1.73 

1.64 

a  =  2.0 

NONE 

-18 

0.29 

2.02 

1.87 

Cauchy 

a  =  3.0 

-31 

-6 

0.42 

1.48 

1.34 

1 

a  =  4.0 

-35 

-6 

0.33 

1.76 

1.50 

.  1+(crX)2  

a  =  5.0 

-30 

-6 

0.28 

2.06 

1.68 

Gaussian 

a  =  2.5 

-42 

-6 

0.51 

1.39 

1.33 

exD(-a2X2/2) a  -  3-0 

-55 

-6 

0.43 

1.64 

1.55 

a  =  3.5 

-69 

-6 

0.37 

1.90 

1.79 

Dolph- 

-50 

0 

0.53 

1.39 

1.33 

Chebyshev 

-60 

0 

0.48 

1.51 

1.44 

-70 

0 

0.45 

1.62 

1.55 

-80 

0 

0.42 

1.73 

1.65 

Kaiser- 

-46 

-6 

0.49 

1.50 

1.43 

Bessel 

-57 

-6 

0.44 

1.65 

1.57 

-69 

-6 

0.40 

1.80 

1.71 

-82 

-6 

0.37 

1.93 

1.83 
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RESOLUTION 


The  reciprocal  relationship  between  aperture  and  beam  width  can  easily  be 
seen  by  examination  of  the  results  from  the  rectangular  window. 

F  (k,o)  )  =  Lsin(kL/2)/kL/2  *  Lsin(7rLsin0/X ) /irLsinG /X  (15) 

&p  o 

The  radiation  pattern  (here  we  use  the  principle  of  reciprocity)  is  simply 

Frt  (k,ai  )  Fftf(kjO)  )  and  represents  the  Farfield  radiation  pattern  sometimes 
fcp  o  £p  o 

referred  to  as  the  Fraunhofer  region  as  opposed  to  Fresnel  and  nearfield 
regions . 


F  (k,U)  ) 
J£p  o 


FI  (k,u>  ) 
£p  o 


L2sin2(TTLsin9/X)/  (ffLsin© /X ) 2 


(16) 


The  result  is  as  expected;  that  for  a  perfectly  coherent  incident  field  an 
increase  in  antenna  length  results  in  additional  power  output  and  a  reduction 
in  beam  width.  These  effects  are  illustrated  in  Figure  3  for  both  space  L  and 
time  T. 


BWFN  =  M  =  2/L/A 

L/A 


(a) 


HPBW  =  50  8  =  0.886 /L/A 

L/A 

(b)  •  • 

1  2 
SOURCES 


Figure  3.  (a)  Power  pattern  and  (b)  power  patterns  for  two  identical  point 

sources  separated  by  the  Rayleigh  angle  as  observed  individually 
(dashed)  and  together  (solid). 
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The  resolution  relationship  can  be  explored  by  a  variety  of  criteria. 

19 

Lord  Rayleigh* s  Criterion  (1897)  was  the  resolvability  of  two  sources  when 
the  maximum  of  the  pattern  for  one  coincided  with  the  first  minimum  of  the 
other.  This  corresponds  to  one  half  the  beam  width  between  first  nulls  (BWFN) 


Rayleigh  Resolution  Criteria  R  88  BWFN/2 


(17) 


The  first  null  occurs  when: 


nLsinO/A  =  it  ;  Lsin0/A  =  1 


(18) 


The  beamwidth  between  the  first  nulls  would  then  be: 


A6fn  =  2A/L 


(19) 


The  alternative  description  also  shown  in  Figure  3  is  to  choose  the  half 
power  (hp)  points  which  can  be  seen  to  be: 


A0,  =  0.886  A/L;  R  =  A/L  =  AG^/2  =  1.12  A0,_ 

hp  FN  hp 


(20) 


This  relationship  holds  for  the  broadside  case  of  the  spatial  array  and 
can  be  generalized  for  the  off  broadside  and  end  fire  cases  as  follows: 


Near  Broadside 
however 
Near  End  Fire 


A0-L/A  =  0.886  sec0 


A0(L/A)1/2  =  2(0.886)1/2 


(21) 


(22) 


this  is  simply  a  consequence  the  properties  of  the  trigonometric  relations. 

The  relations  indicated  here  are  applicable  to  other  window  functions  such  as 

14 

those  shown  in  figure  1  and  can  be  found  in  Harris  as  shown  in  table  1.  The 
major  conclusion  drawn  here  is  that  in  general 


A0-L/A  «  C£  =  1 


(23) 
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Af-At  «  Ct  «  1 


(24) 


represents  the  resolution  properties  of  space-time  processors. 

These  results  are  based  solely  on  the  deterministic  response  of  the  line 
array.  However,  as  was  mentioned  before  in  connection  with  equation  11,  one 
must  require  that  the  resolution  of  the  system  be  less  than  or  equal  to  the 
angular  spread  in  the  incident  field.  Equation  11  shows  that  the  measured  out¬ 
put  is  the  convolution  of  the  system  response  and  the  field  frequency-wavenum¬ 
ber  transform.  In  an  analogous  fashion,  one  can  relate  the  output  of  the  mea¬ 
surement  system  as  the  convolution  of  the  array  response  f  (k,aj)  and  the 

x,p 

Fourier  transform  of  the  acoustic  spatial  correlation  function 


S  (ksin0,u))  =  fn  (ksin0,cj)  *  F(R  (r,*r)). 
P  P  “ 


(25) 


The  acoustic  spatial  correlation  function,  when  we  have  suppressed  the  temporal 
variables,  can  also  be  defined  as 


Rp(ris  jr 2 )  =  (pOr^)  p  (r2>} 


(26) 


where  _r^  and  _r^  are  the  spatial  location  points,  p  is  the  complex  acoustic 
pressure,  t  indicates  conjugation,  and  the  braces  indicate  ensemble  averaging. 
The  relationship  between  R^  and  the  receiving  array  beam  power  output  s  is 
shown  in  equation  27  for  the  case  of  standard  beamforming  in  an  acoustic  field 
that  is  homogeneous  across  the  receiver: 

S(ksin0)  =  / f  (ksin0  -  ksin0*)  {  J R  exp( -iksinG ' r)dr}  d(ksin0’) 

^P  •  P 


=  f0  *  F(R  ) 
Ip  P 


(27) 


In  this  case  the  array  output  as  a  function  of  horizontal  angle  0  is  sim¬ 
ply  the  convolution  of  the  array  beam  pattern  f  with  the  angular  distribution 

£p 
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of  incoming  acoustic  signals,  which  is  given  by  the  spatial  Fourier  transform 
F  of  the  acoustic  correlation  function.  The  symbol  k.  is  the  acoustic  wavenum¬ 
ber  at  the  receiver,  and  the  Fourier  transfom  is  taken  with  respect  to  the  spa¬ 
tial  coordinate  £. 

In  long  range  propagation  cases,  it  is  expected  that  the  initially  smooth 
curve  wave  fronts,  in  the  ray  theory  sense,  will  have  become  irregular  and  spa¬ 
tially  varying  descriptions  of  forward  scattered  energy  bundies.  These  "arriv¬ 
als"  are  resolvable  in  order  but,  when  displayed  after  suitable  steering 
delays,  have  a  phase  and  amplitude  variability  as  well  as  extended  durations 
greater  than  the  transmitted  pulse.  These  distorted  pulses  decorrelate  in 
space--a  phenomenon  equivalent  to  an  angular  redistribution  of  energy. 

The  effects  of  disparities  between  the  array  properties  and  the  incoming 
signal  angular  distribution  are  immediately  apparent  from  equation  27.  A  short 
array  with  a  beam  pattern  broader  than  the  signal  angular  distribution  is 
accepting  unnecessary  noise  on  the  main  lobe,  whereas  a  long  array  with  a  beam 
pattern  narrower  than  the  signal  angular  spread  is  rejecting  part  of  the  power 
in  the  signal.  Thus,  the  design  of  optimal  array  characteristics  is  facili¬ 
tated  by  a  prior  knowledge  of  the  acoustic  correlation  function. 

The  angular  spread  of  arriving  acoustic  waves  during  a  given  observation 
period  can  also  be  due  to  the  arrival  of  multiple  rays  at  different  angles. 

That  is  to  say,  off  broadside  arrivals  from  the  same  azimuthal  angle  but  with 
different  vertical  arrival  angles  can  have  conical  angles  such  that  the  differ¬ 
ence  between  the  conical  angles  from  multiple  rays  is  greater  than  the  half 
power  beamwidth  of  the  array.  Thus  one  will  observe  a  bifurcation  or  the 
appearance  of  acoustic  intensity  on  more  than  one  beam.  The  angle  of  arrival 
for  a  given  path  or  group  of  paths  depend  on  the  relative  orientation  and 
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motion  between  the  source  and  receiver  during  the  observation  period.  Conse¬ 
quently,  motion  both  radially  and  with  respect  to  relative  azimuth  can  have  a 
smearing  effect  on  the  beam  response.  This  phenomena  of  coupling  between  the 
angular  spread  and  the  relative  source-receiver  motion  is  treated  later  in  the 
paper . 

An  analogous  situation  exists  for  the  frequency  domain.  Spread  in  the 
frequency  domain  can  be  caused  by  random  media  scatter,  doppler  spread  due  to 
multipath  projections  of  radial  velocity  and  accelerations.  The  doppler  spread 
due  to  the  relative  motion  source  receiver  is  normally  the  dominant  effect  and 
the  remaining  causes  of  frequency  spread  are  small. 

12 

STATISTICAL  CHARACTERISTICS 

The  determination  of  the  response  of  an  array  to  a  signal  with  variability 

in  phase  is  of  major  importance  to  our  discussions.  The  following  treatment  of 

the  response  of  a  continuous  line  array  to  a  random  signal  closely  follows  the 

12 

derivation  of  Shifrin  and  is  included  here  because  the  results  do  not  appear 
to  be  used  extensively  in  acoustic  antenna  applications  despite  some  previous 
results  worthy  of  note.^  ^  Starting  with  the  transform  of  a  uniform  line 
array  of  length  L  illuminated  by  a  harmonic  wave  p(y)  where  the  time  dependence 
is  suppressed,  one  has 


(28) 


for  convenience,  we  let 


ifj  =  TrLsinG/A  and  x  =  2y/L 


f(^)  -  f  -2/L  =  P(x)exp(itJ;x)dx 


(29) 
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The  pressure  field  p(x)  is  assumed  to  have  a  phase  <j>(x)  which  is  a  zero 
mean  random  variable  with  a  Gaussian  correlation  function. 

P(x)  =  P  exp(i<j)x),  <f>(x)  =  0,  F(x)  =  a  2  =  a  (30) 

o  9 


4>(x,)  <J)(x  ) 

R  =  —7 - x  / - r  =  exp{-(x  -  x  )2/C  2},  C  =  2p  /L  (31) 

<j)  a  (x  )  a  (x  )  1  2  <f>  <{> 


upon  substitution  one  has 

f(ifO  -  PQexp{  i<f)  (x)  4*  ii|;x}dx  (32) 

The  quantity  f(i|0  is  the  spatial  wavenumber  transform  of  the  field.  One 
can  compute  the  mean  wavenumber  transform  and  the  power  wavenumber  spectra 
(which  is  completely  analogous  with  the  farf ield-power-radiation  pattern  of  an 
antenna  with  variable  current  drive  within  the  aperture).  These  computations 

provide  the  basis  for  the  determination  of  the  role  the  variance  (a  2  =  a)  and 

<P 

the  correlation  length  (C,  25  2p/L)  have  on  the  measured  antenna  characteris- 

<P 

tics.  Since  the  process  is  a  stationary  Gaussian  process,  the  averaging  will 
be  ensemble  averaging.  The  mean  field  is 


f(^)|  =  Poexp{i<f>(x)  +  ii|;x}dx  =  P^exp{ 


-a , 2 


+  ii|;x}dx 


=  P  exp 
o 


(33) 


The  important  determination  of  the  quantity  exp(i<f>(x))  is  from  the  charac" 
teristic  function  for  a  Gaussian  distribution  and  can  be  found  in  references  2, 
3,  8  and  12.  Similarly  one  has  for  the  farfield  pattern,  again  using 
reciprocity, 


f  (^  )  |  2  =  Po(x)Po(x^)exp{i((J)(x)-(f>(x^))}exp(i^(x-x^))dxdx^ 


(34) 
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This  integral  can  be  simplified  for  both  the  Gaussian  and  exponentially 
distributed  random  variables.  This  paper  will  simply  outline  the  approach  for 
the  Gaussian  case  as  presented  by  Shifrin. 

|f(i(0|2  -  Pq2  exp{aR^  -  a  +  ii|>(x  “  x^)}dxdx^  (35) 

where 

a  =  a.2,  R.  =  exp{-(x  -  xn)2/C  2},  C  =  C,/i/m  =  2p/i/m L  (36) 

<p  (p  I  <p  m  <p 

00  m 

exp  (a  R  )  =  1  +  l  ^-r  exp{  -  (x  -  x  )2/C  2} 

<J>  -  m :  l  m 

m=l 


then 


f(^)|2  “  exp(-a){P2  JJ  exp(  iip  (x-x^ )  )dxdx^  + 


00  m  , ,  /  -(x-x1  )2  \ 

£  P0  ^-1  exp(  - 5 -  +  i*(x“xl))  dxdxl 

m—  1  ' 


m 


(37) 


finally  with  =  1/2  one  has 


f(’l’)  I  2  “  exp(-a){ S1?2  ~  +  1M  l  ~T  KC  ,1^)} 


m=l 


m!  m 


(38) 


The  mean  power  pattern  has  thus  been  shown  to  depend  on  the  pattern  in  the 

absence  of  phase  randomness  and  a  distortion  in  the  pattern  due  to  the  combined 

influence  of  the  phase  variance  a  and  the  normalized  correlation  length,  C  . 

m 

The  power  in  the  direction  of  the  principal  maxima  is  also  seen  to  be  decreased 

by  the  factor  exp(-a).  The  function  l(C  ,ip)  can  be  evaluated  in  closed  form 

m 

for  special  cases.  A  similar  function  can  be  derived  when  the  exponential  form 
of  the  correlation  function  is  used. 
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The  figure  4  from  Shifrin  shows  the  mean  power  radiation  pattern  |f(i^)|2 
versus  the  angle  ip  for  various  values  of  the  phase  variance  a  and  the  normal¬ 
ized  correlation  length  C.  These  calculations  by  Shifrin  show  some  interesting 
general  characteristics.  The  randomicity  in  phase  quantified  by  the  variance 
shows  that  the  effect  of  an  increase  causes  the  principal  maxima  to  decrease 
(see  a,  b,  c,  d  and  e).  In  addition,  the  phase  variance  results  in  a  smoothing 
or  smearing  of  the  radiation  pattern,  that  is,  the  nulls  are  filled  or  blurred. 
As  the  phase  variance  increases,  the  pattern  loses  detailed  structure  and 
becomes  a  monatonically  decreasing  one  (see  a  and  b  for  a  =  3).  The  effect  of 
increasing  the  phase  correlation  length  is  shown  in  figure  4e.  As  C  increases, 
the  pattern  approaches  the  ideal  pattern,  the  power  in  the  direction  of  the 
principal  maxima  increases  (see  example  e  and  f).  In  summary,  the  increase  in 
phase  variance  causes  a  reduction  in  power  for  the  principal  maxima  and  a  blur¬ 
ring  of  the  pattern.  An  increase  in  the  phase  correlation  length  causes  the 
power  for  the  principal  maxima  to  increase  and  the  pattern  to  approach  the 
ideal  case.  These  were  numerical  computations  performed  by  Shifrin.  The 
expression  can  be  simplified  to  approximate  forms  as  discussed  next. 

When  the  correlation  distance  of  the  phase  fluctuations  is  sufficiently 
great  to  make  the  variation  in  phase  over  a  wavelength  along  the  array  small 
with  respect  to  unity  then  the  mean  directional  gain  can  be  shown  to  be 

00  m 

D  =  P  exp(-a){ 1  +  1/4  7  —  l(C  ,0)}  (39) 

o  m!  m 

m=l 

It  is  instructive  to  consider  the  behavior  of  D  when  the  length  of  the 
array  (L)  is  much  smaller  than  the  correlation  distance  (2p  >>  L,C  >1).  We 
can  further  divide  this  case  into  two  regions.  Region  I  where  C//7  =  > 

1  and  region  II  where  C//a  <  1. 
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iwi 2  (dB) 


|WT|Z  (dB) 


(a) 


|W  (dB) 


(b)  (d)  (f) 


Figure  4.  The  mean  radiation  pattern  |f(i|<)|^  for  various  values  of  the  phase  variance  (a  =  0^) 

and  the  phase  correlation  length  (C  =  2p/VmL)  versus  \p  =  irLsinQ/A.  (a-d)  show  the  mean 
radiation  pattern  of  a  linear  system  for  various  a  with  (a)  C  =  0.25,  (b),  C  =  0.5, 

(c)  C  =  2.0,  (d)  C  =  6.0.  For  all  values  of  C  the  effect  of  a  is  to  decrease  the  intensity 
at  the  principal  maximum  and  to  smooth  the  pattern  by  filling  the  nulls.  (e)  shows  for  a 
constant  a,  the  efffect  of  increase  in  phase  correlation  length  on  the  character  of  the 
pattern  while  (f)  shows  for  constant  correlation  length  the  effect  of  phase  variance 
increase  on  the  principal  maximum. 
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In  the  first  region,  figure  5,  the  mean  directional  gain  can  be  shown  to 
be 


C/'/a  >  1,  D  =  D  I(C  ,0)  =  D  (1  -  2a/3C2)  =  D  (40) 

o  a  o  o 

which  is  equal  to  the  ideal  gain  of  the  line  array.  In  the  second  region,  the 
directional  gain  can  be  approximated  by 

C/Ja  <  1,  D  «  D  I(C  ,0)/4  *  D  /2tF  C/4/T  =  2»^7a(p/X)  (41) 

O  CL  O 

This  expression  shows  the  mean  directional  gain  to  be  independent  of  array 
lengths,  i.e.,  saturated. 


Figure  5.  Mean  directive  gain  versus  length  (L/A). 

One  may  interpret  these  results  by  considering  the  relationship  between 
the  phase  variance  a  =  A(J)2  and  the  variance  of  the  arrival  angle  0  of  the 
acoustic  wave  with  respect  to  the  array.  Since  p  >  L,  the  mean  phase  variation 
across  the  array  may  be  taken  as  A<f>  =  kL0 .  For  this  case  of  p  >  L  we  have  for 
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the  Gaussian  variable  A<f>2  =  2a  L2/p2.  Equating  these  quantities  leads  to  the 
result  for  the  mean  square  variation  in  the  arrival  angle  0  to  be 

02  =  8(a  L/p  )2/k2L2  =  8(i/a/c)/k2L2  (42) 

Then  this  estimate  of  the  angular  spread  in  the  acoustic  field  can  be 
compared  to  the  half-power  width  of  the  antenna  pattern. 


0,  =  20,  =  A/L;  (C/r£)2  =  20,  /(tt0)2 

hp  hp  hp 


(43) 


When  0^  <  92  then  the  spread  in  the  angle  of  the  incident  energy  is 

greater  than  the  half-power  beam  width  and  <  1 .  In  this  instance,  our 

antenna  exhibits  saturation.  On  the  other  hand,  when  ©  >  02,  the  beam  width 

hp 

of  the  antenna  is  greater  than  the  incident  spread  of  energy  and  an  increase  in 
length  results  in  an  increase  in  mean  directional  gain. 


The  final  case  is  when  the  length  of  the  antenna  is  much  greater  than  the 

correlation  length  (L  >  p,  c  <  1).  In  this  instance,  equation  39  for  the  mean 

directive  gain  becomes  (I(C  ,0)  •+  27rc/v/m) 

m 

nr  °°  m 

D  =  D  exp(-a){l  +  — -  V  7-7=}  =  D  exp(-a)  (1  +  CN(a))  (44) 

o  4  u  -  m !  /m  o 

m=l 


when  CN(a)  >  1,  this  term  dominates  the  gain  and  since  (D  =  D  exp(-a)CN(a)  = 

o 

4pexp(a)N(a)/4  =  D  is  again  independent  of  array  length,  the  array  is  still 
saturated.  In  the  case  of  radar  waves  in  the  atmosphere,  Shifrin  maintains 
that  further  increases  in  the  length  of  the  system  produces  a  point  where  CN(a) 
<  1  and  then  the  mean  directive  gain  becomes  D  =  D^e  a.  That  is  although 
reduced  by  the  constant  factor  e  a  the  mean  directive  gain  increases  in  propor¬ 
tion  to  the  ideal  gain,  -  2L/X .  This  region  III  which  is  dominated  by  the 
correlation  length  is  shown  on  figure  5.  The  existence  of  a  region  III  for  the 
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non-plane  multipath  underwater  acoustic  propagation  case  is  still  questionable. 
For  wavelengths  greater  than  1.5  m,  practical  measurements  in  the  ocean  will  be 
limited  to  log(L/A)  _<  4.  The  region  II  can  be  bounded  by  letting  the  approxi¬ 
mate  expression  for  CN(ot)  =  1.  When  this  is  done,  one  arrives  at  an  upper 
bound  of  region  II 


00 


and  the  lower  bound  of  region  II 


(45a) 


L 


Z 


2p  /  VcT 


(45b) 


Thus  the  array  behavior  can  be  characterized  provided  one  has  a  knowledge 
of  a  and  p . 


In  summary,  the  presence  of  phase  randomicity  was  shown  to  have  marked 
effects  on  the  mean  radiation  pattern  and  mean  directive  gain.  The  variance 
was  found  to  cause  a  blurring  of  the  pattern,  i.e.,  the  nulls  are  blurred  and 
the  principal  maxima  decreases.  As  the  variance  increases  further  the  pattern 
was  shown  to  become  a  monatonically  decreasing  function.  As  the  correlation 
distance  was  increased  the  pattern  was  found  to  approach  the  characteristic  of 
a  system  without  phase  fluctuations.  The  mean  directional  gain  was  shown  to 
exhibit  three  types  of  behavior:  the  firsts  a  region  in  which  the  gain  in¬ 
creases  proportional  to  2L/A;  the  second,  a  region  of  gain  saturation  (gain 
independent  of  length);  and  the  third,  in  which  the  gain  increase  by  2L/A 
attenuated  by  exp(-a).  (This  is  questionable  for  long  range  underwater  acous¬ 
tic  propagation.)  These  results  were  derived  for  a  Gaussian  correlation  func¬ 
tion  but  could  also  be  developed  for  an  exponential  distribution. 
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In  closing  this  section,  the  Gaussian  correlation  function  is  related  to 
the  measure  of  spatial  coherence.  Recalling  equations  30,  31,  and  35,  one  can 
indicate  this  relationship  by  noting 

R^  =  exp{  -  (x  -  Xj)2/C^2}  -  1  -  |  x  -  x^ | 2/C2  (46) 

substitution  of  this  quantity  into  equation  35  results  in 

|f(^)|2  =  JJ+J  P2exp{-(x  -  x^)2/C2/2X}exp(ii/;(x  -  x^))  dxdx^  (47) 

one  recognizes  P2exp( -(x-x2 )/L2)  as  a  Gaussian  form  of  the  spatial  coherence 
O  i  H 

function. 
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COHERENCE 


The  statistical  response  of  an  array  has  been  developed  for  a  signal  with 
a  random  component  of  phase  which  is  Gaussian,  Due  to  general  characteristics 
of  the  Gaussian  process,  it  can  be  shown  that  a  random  phase  variable  with  a 
Gaussian  spatial  correlation  function  also  produces  a  spatial  coherence  func¬ 
tion  which  was  found  to  have  a  Gaussian  form.  Although  not  shown  in  this 
paper,  the  random  phase  variable  with  an  exponential  correlation  function  pro¬ 
duces  an  exponential  spatial  coherence  function.  The  results  have  primarily 
dealt  with  the  spatial  dimension  but  also  apply  to  the  temporal  scale. 

The  analogous  nature  of  the  treatment  of  coherence  in  space  and  time  is 
restated  in  Figure  6.  The  top  set  of  expressions  the  spatial  coherence  whereas 
the  bottom  sets  describe  the  temporal  coherence. 


R 


p 


SPATIAL  COHERENCE 

(X)  P*  (X  +  RJ}>  Mutual  coherence  function 


Fs  W 

S(k  Sin  6)  =  BA  *  Fs  {Rp} 


Array  beam  signal  output 


Incident  signal  angular  distribution 


TEMPORAL  COHERENCE 


Q 


Autocovariance  function  or  fluctuation 
time  series 


Ft  {q} 

W  *  Ft  {Q} 


Processor  signal  output 


Incident  signal  frequency  spectrum 


Figure  6.  Spatial  and  temporal  coherence  relationships. 
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The  top  equation  defines  the  spatial  covariance  of  the  acoustic  field,  R  , 

P 

to  be  the  average  of  the  product  of  the  complex  pressure  at  one  point  in  the 
space  domain  with  the  complex  conjugate  of  the  pressure  at  a  second  point  in 
space.  In  a  number  of  theories  for  wave  propagation  in  a  random  medium,  this 
expression  is  called  the  mutual  coherence  function.  The  spatial  Fourier  trans¬ 
form  of  the  mutual  coherence  function  gives  the  incident  signal  angular  distri¬ 
bution.  It  was  shown  for  the  case  of  standard  beamforming  and  an  acoustic 
field  that  is  homogeneous  across  the  aperture  that  the  array  beam  signal  out¬ 
put,  S(oo,k),  is  the  convolution  between  the  array  beam  pattern  and  the  incident 
signal  angular  distribution. 

The  bottom  set  of  expressions  treats  temporal  coherence.  The  first  equa¬ 
tion  defines  the  temporal  autocovariance  function,  Q,  as  the  average  of  the 
product  of  the  complex  pressure  at  one  point  in  time  with  the  complex  conjugate 
of  the  pressure  at  another  point  in  time.  This  is  sometimes  referred  to  as  the 
autocovariance  function  for  the  fluctuation  time  series.  The  temporal  Fourier 
transform  of  Q  gives  the  incident  signal  frequency  spectrum.  It  can  be  shown 
that  the  processor  signal  output  is  the  convolution  of  the  processor  spectral 
window  with  the  incident  signal  frequency  spectrum. 

SPATIAL  COHERENCE 

The  propagation  of  waves  though  a  random  media  has  been  treated  by  Cher- 
20  21 

nov  and  Tatarski  .  The  results  of  Chernov* s  work  and  the  experimental  mea- 

12 

surements  discussed  by  Shifrin  were  the  basis  for  Shifrin's  choice  to  use 
exponential  and  Gaussian  phase  correlation  functions  to  determine  the  statisti¬ 
cal  characteristics  of  antennas.  The  Gaussian  form  was  chosen  to  illustrate 
the  general  characteristics  of  line  array  when  the  incident  wave  had  a  normally 
distributed  phase  variable  with  a  Gaussian  correlation  function.  The  exponen- 
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tial  correlation  function  was  not  chosen,  even  though  in  several  instances  this 
function  compares  favorably  with  experimental  data  on  electromagnetic  propaga¬ 
tion  in  the  atmosphere,  since  an  exponential  correlation  function  represents  a 
discontinuous  physical  process.  The  results  of  this  treatment  starting  with 
the  Gaussian  phase  correlation  function  were  shown  to  produce  an  equivalent 
result  with  a  Gaussian  spatial  coherence  function  (see  eq.  47). 


The  ocean  acoustics  wave  propagation  problem  is  thought  to  be  more  complex 
than  atmospheric  propagation  due  to  the  spatial  properties  which  are  anisotrop¬ 
ic  and  nonhomogeneous .  A  thorough  discussion  of  this  subject  is  beyond  the 

22-27 

scope  of  this  paper  and  we  refer  the  reader  to  some  very  excellent  reviews. 

(This  subject  has  been  extensively  studied  and  comprehensively  reviewed  by 

22  23  25  26 

Flatte.  5  Desaubies  has  prepared  a  readable  overview  while  McCoy  and 

27 

Guthrie  have  reviewed  the  work  performed  at  the  Naval  Research  Laboratory.) 


We  have  shown  that  the  output  of  a  line  array  can  be  written  as 


S(ip)  =  £  (ty-ty  )  *  F(R  ) 

£p  o  p 


(50) 


The  spatial  coherence  function  as  shown  in  figure  6  can  be  written  as 


R  =  <P(x)  P+(x  +  r ) >  =  <P2>  exp(-(y/L  )n)  n  =  1,  1.5,  2  (51) 

P  ~  ~  “  h 

Where  <P2>  is  the  mean  square  pressure  and  is  proportional  to  the  acoustic 
intensity  and  where  n  =  1  is  the  exponential  form,  n  =  1.5  the  Beran-McCoy 
form,  and  n  =  2  the  Gaussian  form.  Since  our  goal  is  to  understand  line  array 
performance  in  the  linear  region  of  figure  5  to  the  beginning  of  the  gain  satu¬ 
ration  region,  these  types  of  coherence  functions  can  be  considered  valid. 

22  23  24 

Flatte  5  and  Dashen  relate  the  coherence  to  a  phase  structure  function. 
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R  =  <P2>  exp(-D(l,2)/2)  (52) 

P 

and  in  the  region  of  interest  here,  the  geometrical  acoustics  and  partially 
saturated  ($A  ^  1)  regions,  they  give  the  structure  function  at  small  separa¬ 
tions  to  be 

D(l,2)  =  $2{(T/t)2  +  (y/LR)2  +  (Z/Lv)2(Log(Z/Lv)  +  ry)>  (53) 

The  characteristic  time  t  =  1.6  hours;  the  characteristic  vertical  coher¬ 
ence  length,  L^,  is  300  m;  and  the  characteristic  horizontal  coherence  length, 

L  ,  is  3. 7-6. 4  km.  The  structure  function  given  here  clearly  corresponds  to  a 
H 

coherence  function  with  n  =  2,  the  Gaussian  form. 

28 

The  exponential  functional  form  of  the  coherence  was  employed  by  Cox 
since  it  enabled  the  development  of  a  closed  form  solution  to  the  array  gain 
problem.  Cox  used  the  function  only  for  small  (A/2)  separations  and  the  result 
is  useful.  However,  as  Shifrin  has  pointed  out,  this  functional  form  has  a 
finite  derivative  at  zero  separations  and  corresponds  to  a  discontinuous  physi¬ 
cal  process,  not  representative  of  wave  propagation. 

The  work  presented  in  this  paper  is  based  on  the  n  =  1.5  functional  form 

29  30 

after  the  work  of  Beran  and  McCoy.  5 

Theoretical  predictions  for  the  acoustic  correlation  function  can  be 
obtained  by  solution  of  a  Bethe-Salpeter  equation  that  includes  multiple  scat¬ 
tering  and  the  effects  of  random  oceanic  fluctuations: 

R  ( X-  ,  X  )  =  R  (X,  XrJ 
p  —1  —2  p  —  — 0 

ro 

+  ff  G(x2,P)  k2(p_)  G(p_,p_!)  k2(p_')  Re(p_,P_’)  R^(x,P_)dpdp  ’  (54) 
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+  JJ  G(xlfr)  k2(r)  G(r,r')  k2(r')  R^r.x')  R  (r' ,x2)drdr ' 

+  JJ  G(x^,r)  k2(r)  GCx^.fO  k2(fO  R^Cr^fO  Rp(_r  ,p_)drdp 

31 

The  Bethe-Salpeter  equation  (Frisch  (1968))  given  in  equation  54  is 

appropriate  for  a  weak  environmental- fluctuation  field  and  is  truncated  at  the 

two-point  moment  of  the  environment  fluctuation.  This  integral  formulation 

indicates  that  the  acoustic  correlation  function  at  the  receiver  is  a  function 

of  the  following:  R  ,  the  correlation  function  of  the  solution  of  the  reduced 

Po 

deterministic  wave  equation;  G,  the  Green* s  function  for  the  reduced  determin¬ 
istic  wave  equation;  k,  the  deterministic  wavenumber;  R^  ,  the  two-point  corre¬ 
lation  function  for  the  random  portion  of  the  wavenumber;  and  R  ,  the  acoustic 

P 

correlation  function  in  the  space  between  the  source  and  receiver. 

Four  of  the  usual  assumptions  made  to  provide  tractability  to  equation  54 
are  that  (1)  G  is  the  infinite  space  Green's  function;  (2)  the  variation  of  the 
deterministic  wavenumber  is  parametric  with  range;  (3)  the  acoustic  field  is 
statistically  homogeneous  over  distances  that  are  large  compared  with  the 
receiver  aperture;  and  (4)  the  environmental  field  is  statistically  homogeneous 
and  isotropic.  While  none  of  these  assumptions  is  strictly  valid  for  long- 
range  propagation,  the  first  three  assumptions  may  be  considered  approximately 
valid  in  a  perturbation  sense.  However,  the  inherent  anisotropy  and  statisti¬ 
cal  inhomogeneity  of  the  index  of  refraction  in  the  vertical  direction  cause 
the  fourth  assumption  to  be  extremely  questionable. 

29-31 

Beran  and  McCoy  derived  a  solution  in  which  the  basic  anisotropic 

nature  of  the  two-point  environmental  statistic  was  retained  together  with  the 
approximate  form  of  the  power  spectrum  of  random  temperature  fluctuations  in 
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the  horizontal  plane.  One  of  the  results  of  the  Beran-McCoy  formulation  is 
given  as  follows 

R  (r,f,y,z)  =  I ( z )  exp(-E  f^2ry^2)  =  j(z)  exp(-E  k^2ry^^)  (55) 

p  r  k 

In  this  equation,  the  acoustic  correlation  function  along  a  horizontal 
line  transverse  to  the  direction  of  propagation  is  an  exponential  function 
whose  argument  is  proportional  to  E ,  an  environmental  parameter  appropriate  for 
the  propagation  path;  the  frequency  f  to  the  5/2  power;  the  range  R  and  the  3/2 
power  of  Y  the  transverse  separation  distance  of  the  two  correlation  points. 

The  symbol  I  represents  the  intensity  at  a  single  point  receiver  with  the  same 
range  and  depth  Z. 


The  components  of  the  environmental  parameter  E  for  the  finite  source  case 
are  found  by  Beran  and  McCoy  to  be 


E  =  1.7C1/C 
f  o 


r)  c 

g^)2A  L  (21T/C  )5/2 
d  T  I  ym  o 


(56) 


where  is  the  nominal  speed  of  sound,  SC^/ST  is  the  partial  derivative  of  the 
sound  speed  with  respect  to  the  temperature,  A^,  is  the  coefficient  of  the  sin¬ 
gle  term  power-law  spectrum  representation  of  the  random  horizontal  temperature 
variations  (i.e.,  the  nominal  strength  of  the  random  temperature  field),  and 
L  is  the  correlation  length  of  the  random  temperature  fluctuations  in  the 
vertical  direction.  Thus  the  spatial  coherence  problem  can  be  considered  the 
experimental  determination  of  the  functional  form  of  the  coherence  and  the  mea¬ 
surement  of  the  coherence  length.  These  specify  the  performance  of  the  line 
array . 


*Note:  If  the  expression  for  R 

P 

instead  of  frequency  f  we  have  E 


(eq.  55)  is  written  in  terms  of  wave  number  k 

,  =  (C/2tO"5/2  E,  . 
f  k 
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SPATIAL  COHERENCE  AND  THE  LINE  ARRAY 

Since  the  response  of  the  array  to  the  sound  field  is  the  convolution  of 
the  beam  response  of  the  array  and  the  wavenumber  spectra  of  the  field,  we  can 
establish  a  resolution  criteria  by  matching  the  angular  spread  in  the  incident 
field  to  the  3  dB  width  of  the  beam  response.  The  spatial  coherence  function 
is  a  measure  of  this  angular  spread  and  can  be  evaluated  as 

S(ksin0)  =  j  m  R(y)  exp{ iksin0y}dy  (57) 

=  J_oo  exp( - (y/L^)n)  exp(iksin0y)dy 

This  integral  can  be  evaluated  in  closed  form  for  n  =  1,  2  and  numerically- 
for  n  =  1.5.  The  results  of  this  analysis  when  matched  to  the  3  dB  width  yield 
the  following: 


n 

BW  (RAD) 

BW  (deg) 

L  /L. 
a 

1 

0.318  A/L 
0.457  H 

18.2°  A/L 
26.2° 

2.72 

1.5 

1.89 

2.0 

0.530 

30.36° 

1.64 

f*P 

0.886 

50.76° 

1 

That  is  the  exponential  form  has  an  angular  spread  of  0.318  A/L  and  when 

H 

matched  with  the  3  dB  beam  width  of  the  line  array  0.886  A/L  yields  an  aper- 

3. 

ture  length  2.72  times  the  horizontal  coherence  length.  Note  that,  given  the 
same  coherence  length,  the  acoustic  field  angular  spread  for  the  Gaussian  (n  = 
2)  form  is  1.67  times  as  large  as  the  acoustic  field  angular  spread  for  the 
exponential  (n  -  1)  form. 
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In  addition  to  matching  the  angular  spread  in  the  incident  field  to  the 
beam  width,  one  can  compute  the  relative  signal  gain  (RSG),  that  is,  the  ratio 
of  actual  signal  gain  to  ideal  signal  gain.  It  can  easily  be  shown  that 


(58) 


a'  H 


This  relationship  shows  that  the  relative  gain  of  the  system  (when  related 


to  an  exponential  form  of  the  spatial  coherence  function)  is  a  function  of  the 


coherence  length.  Thus  the  angular  resolution  (beamwidth)  and  gain  of  the  line 
array  are  determined  by  the  coherence  length  and  the  form  of  the  coherence 
function.  This  integral  is  calculated  for  and  the  results  are  shown  in  figure 
7.  The  relative  signal  gain  is  plotted  versus  the  length  of  the  array  in  units 
of  coherence  lengths.  The  limiting  value  of  n  =  1,  1.5,  2  corresponding  to  the 
matched  beam  width  ratio  of  1.93  for  the  Beran-McCoy  form  is  also  shown.  The 
experimental  determination  of  coherence  length  using  the  functional  form  shown 
in  equation  55  allows  one  to  bound  the  limits  on  array  length  due  to  the  ocean 
volume . 

EXPERIMENTAL  RESULTS 

To  compare  the  theoretical  correlation  function  with  the  acoustic  measure- 
32 

raents  by  Stickler  ,  one  must  ascertain  the  value  of  the  environmental  parame¬ 
ter  E  that  is  appropriate  for  the  experimental  situation.  The  terms  involved 

in  the  computation  of  E  were  obtained  in  the  following  manner.  Because  most  of 

32 

the  acoustic  data  were  based  upon  the  transmission  path  that  intersected  the 
array  at  approximately  a  13°  elevation  angle,  a  ray  tracing  program  was  imple¬ 
mented  to  describe  the  path  (range  versus  depth)  of  the  13°  ray  from  the 
receiver,  using  measured  sound  speed  profiles  along  a  radial  from  the  array  at 
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ARRAY  LENGTH  IN  UNITS  OF  COHERENCE  LENGTH 

Figure  7.  Relative  signal  gain  versus  array  length  in  units  of  coherence 
length  for  different  coherence  forms. 

the  azimuth  of  the  experiment.  The  average  sound  speed  along  the  13°  ray  path 
was  1517.5  m/ s ,  and  cq  was  set  equal  to  this  value.  The  constant  3.6  m/s/°C 
was  the  value  used  for  the  partial  derivative  of  the  sound  speed  with  respect 
to  temperature,  and  this  was  obtained  by  evaluating  the  change  with  respect  to 
temperature  of  Wilson's  sound  speed  equation  at  a  nominal  temperature  of  10°C. 


The  procedure  that  determined  was  as  follows.  First,  the  Brunt-Vaisala 

frequency  N  was  averaged  along  the  13°  ray  path,  and  the  average  value  was 

■3  33 

found  to  be  1.5  x  10  (rad/s).  Then  the  measured  coefficients  to  the  single 

term  power-law  spectra  for  the  random  horizontal  temperature  fluctuations  were 

treated  as  a  function  of  N.  A  linear  extrapolation  of  the  measured  coefficient 
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from  the  nearest  measured  N  to  the  average  N  along  the  13°  ray  path  indicates  a 
value  of  A of  1.5  x  10  ^  °C^/m. 

The  normalized  temperature  spatial  correlation  function  is  equal  to  the 
inverse  Fourier  transform  of  the  temperature  spatial  power  spectrum  normalized 
by  the  total  variance;  and  the  correlation  length  can  be  defined  as  the  spatial 
length  for  which  the  value  of  the  normalized  correlation  function  is  0.5. 

Thus,  the  correlation  length  in  the  vertical  direction  of  the  temperature  fluc¬ 
tuations  was  estimated  by 

0.5  -  fk  F  (k)  cos(ks)  dk/f^  F  (k)  dk  (59) 

J  a  v  J  a  v 

where  F^  is  the  vertical  spatial  power  spectrum  of  the  temperature  fluctua¬ 
tions,  k  is  the  vertical  wavenumber,  s  is  the  estimate  of  L  ,  and  a  =  .0209 

ym 

34 

rad/m  and  b  =  20.94  rad/m.  Millard  has  data  taken  southwest  of  Bermuda  which 

for  the  wavenumber  interval  under  consideration  can  be  approximated  by  F^(k)  = 

ck^,  where  d  is  between  2  and  2.5.  The  value  for  s  is  *  35  m  in  the  case  when 

d  is  2,  s  is  -  28  m  when  d  is  equal  to  2.5.  It  might  also  be  noted  that  the 

35 

theory  of  Garrett  and  Munk  would  predict  d  to  be  equal  to  2.5.  As  a  result 

of  these  considerations,  L  was  given  the  nominal  value  of  30  m.  Insertion  of 

ym 

the  values  of  into  equation  58  yields  the  environmental  parameter  value 

E  =  4.8  x  10"17  (60) 

This  environmental  parameter  was  also  calculated  for  sound  speed  profiles  from 
the  Mediterranean,  Pacific  and  the  Arctic  and  are  summarized  as  below: 


Dm 

E 

L 

c 

Atlantic 

4000 

m 

4.88  x  1017 

1634 

m 

Mediterranean 

2000 

m 

4.88 

1634 

m 

Pacific 

3700 

m 

9.31 

1042 

m 

Arctic 

1600 

m 

13.06 

793 

m 
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Dm  is  the  maximum  depth  over  which  the  parameter  E  was  averaged.  The 
smallest  E  were  for  the  Atlantic  and  Mediterranean  whereas  the  largest  was  for 
the  Arctic.  The  variation  in  E  is  2.7  with  a  coherence  length  variation  of  2. 
That  is,  the  coherence  length  does  not  change  by  orders  of  magnitude  as  would 
be  indicated  without  averaging  the  environment  along  the  propagation  path. 

To  compare  measured  acoustic  results  with  the  theoretical  form  of  the  cor¬ 
relation  function  without  having  a  variable  parameter  E  involved  in  each  com¬ 
parison,  it  is  necessary  to  obtain  an  acoustic  data  base  wherein  the  propaga¬ 
tion  is  along  environmentally  similar  paths.  In  addition,  because  the  Beran- 
McCoy  theoretical  solution  is  strictly  applicable  only  to  single  path  propaga¬ 
tion,  path-resolving  experimental  capability  is  required. 

The  work  of  Stickler  and  Broek  using  a  planer  receiver  in  combination  with 

short  (10-ms  ping)  acoustic  signal  structure  provided  the  required  path 
32 

resolution . 

During  the  acoustic  measurements,  steering  delays  correcting  for  average 
wave  front  curvature  were  constant  for  each  experimental  sequence  that  had  time 
durations  varying  from  8  to  24  hr.  The  steering  delays  compensating  for  source 
motion  were  allowed  to  vary  linearly  with  time  during  each  experimental  data 
sequence . 

The  acoustic  data  resulting  from  the  experiments  were  reported  as  the 
average  of  the  cosine  of  the  phase  difference  between  pairs  of  receivers  and 
are  plotted  versus  receiver  separation.  A  direct  comparison  between  the  acous¬ 
tic  correlation  function  and  the  average  of  the  cosine  of  the  phase  differences 
was  justifiable  because  data  indicated  that  above  100  Hz  the  phase  fluctuations 
account  for  the  predominant  portion  of  the  change  in  the  acoustic  correlation 
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for  single  path  transmission.  It  should  be  noted  that  although  some  of  the 
references  included  data  for  vertical  separation  between  receivers,  only  data 
for  horizontal  separation  between  receivers  (and  transverse  to  the  propagation 
direction)  are  used  in  this  paper. 

The  Stickler  data  are  shown  in  figure  8  and  figure  9  for  the  center  fre¬ 
quency  of  400  Hz  and  ranges  of  137  km  and  259  km.  We  observe  that  the  spread 
in  the  measured  coherence  values  increases  as  the  separation  between  receivers 
increases.  Note,  the  median  value  of  coherence  decreases  as  the  separation 
increases;  and  for  a  fixed  signal  to  noise  ratio,  small  values  of  coherence  are 

difficult  to  measure  and  can  be  corrupted  by  the  noise.  These  effects  are  con- 

36  37  38 

sistent  with  the  analysis  of  Devilbiss  et  al.  and  Carter.  5  The  results 
can  be  interpreted  as:  the  measurement  of  phase  in  a  random  additive  noise 
background  requires  a  high  signal- to-noise  ratio.  Even  at  6  dB  the  uncertainty 
in  the  phase  may  be  governed  by  the  properties  of  the  additive  noise.  Theore¬ 
tical  analysis  of  the  measurement  of  the  coherence  of  a  fluctuating  signal  in 
an  additive  random  noise  background  also  shows  that  the  coherence  is  a  function 
of  signal  to  noise  ratio,  with  large  confidence  bounds.  The  conclusion  drawn 
is  the  estimation  of  signal  phase  and  coherence  properties  in  the  presence  of  a 
noise  field  such  as  the  ocean  ambient  noise  field  is  difficult.  This  conclu¬ 
sion  is  underscored  by  the  spread  in  the  data  obtained  by  Stickler.  The  smooth 
curves  in  figures  8  and  9  are  the  calculations  performed  using  the  environmen¬ 
tal  parameter  E  and  the  McCoy  formulism.  In  figure  10,  the  theoretical  curve 
is  repeated  and  the  curves  for  all  three  exponential  forms  for  the  coherence 
function  are  shown  where  the  correlation  length  was  determined  from  the  esti¬ 
mated  (1/e)  value  observed  in  the  data.  The  calculation  utilizing  the  estimat¬ 
ed  E  value  seems  to  bound  the  data  nicely  especially  if  one  assumes  the  down¬ 
ward  spread  in  the  coherence  measurements  could  be  due  to  the  additive  noise. 


36 


TR  8077 


Figure  8,  Acoustic  correlation  ver¬ 
sus  transverse  receiver  separation: 
experiment  three;  range,  137  km; 
frequency,  400  Hz;  curve  represents 
theory. 


Figure  9.  Acoustic  correlation  ver 
sus  transverse  receiver  separation: 
experiment  one;  259  km;  frequency, 
400  Hz;  curve  represents  theory. 


R/X=  3.65  •  1 04 


Figure  10. 


A  comparison  of  exp(-(y/L  )n)  with  acoustic  correlation  versus 
transverse  separation  at  £00  Hz. 
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Further,  using  the  empirical  estimate  of  coherence  length,  the  n  =  3/2  form 
provided  a  nice  match  to  the  data  bounded  on  the  upward  side  by  the  n  =  2  form 
and  the  lower  side  by  the  n  =  1  form. 


Based  on  this  comparison,  we  show  in  figures  11  and  12  the  comparison  of 

the  calculated  coherence  versus  range  and  coherence  length  versus  range  at  400 

32 

Hz  based  on  the  measurements  of  Stickler  and  Broek.  Note  that  the  coherence 

lengths  here  are  between  350  to  1647  m  compared  to  the  value  of  3.4  to  6.4  km 
24  39 

used  by  Dashen.  Ross  Williams  has  estimated  coherence  lengths  at  various 

ranges  by  visually  smoothing  phase  time  histories.  His  estimates  ranged 

between  331  m  and  1200  m  at  400  Hz  for  ranges  between  172  and  496  km  as  shown 

41 

by  the  circles  in  figure  11.  Finally,  Kennedy  estimates  the  phase  coherence 
length  to  be  183  m  at  800  Hz  and  a  range  of  46  km.  The  Williams  values  are 
seen  to  be  comparable  to  the  St ickler-Broek  values.  Using  the  functional 
dependence  on  frequency  to  5/2  power  and  range,  we  estimate  that  the  Kennedy 
results  would  correspond  to  a  length  of  approximately  353  m,  the  lower  end  of 
our  range.  Nevertheless,  the  spread  in  the  data  stresses  our  conclusion  that 
the  pairwise  coherence  is  difficult  to  measure  and  is  strongly  influenced  by 
the  signal- to-noise  ratio. 


Fig.  11.  Acoustic  correlation  versus 
range;  fixed  receiver  sepa¬ 
ration,  610  m;  frequency, 
400  Hz. 


Fig.  12.  C  versus  range:  frequency, 
400  Hz;  curve  represents 
theory . 
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To  underscore  this  conclusion,  magnitude  squared  coherence  data  is  pre¬ 
sented  in  figure  13  for  two  frequencies,  175  Hz  from  a  moored  source  and  173  Hz 
from  a  towed  source.  The  MSC  curve  for  the  173  Hz  tone  shows  a  pronounced  mul¬ 
tipath  interference  pattern  comparable  to  the  signal- to-noise  (S/N)  variation. 
The  results  for  the  175  Hz  tone  at  high  S/N  (>  10  dB)  replicates  the  S/N  varia¬ 
tions  with  increasing  separations.  These  measurements  were  performed  with  a 
streamer  behind  a  tow  ship.  Hence,  as  one  approaches  the  ship,  the  interfering 
tow-ship  radiated  noise  increases,  the  S/N  decreases  (<  10  dB),  and  the  MSC 
shows  wide  variations  due  to  the  increase  in  the  percentage  of  interfering 
noise.  This  result  underscores  the  conclusions  drawn  that  pairwise  coherence 
is  difficult  to  measure,  is  strong  function  of  signal- to-noise  ratio,  and  its 
bounds  are  large.  The  estimation  of  coherence  must  also  be  performed  in  a  time 
that  is  short  compared  to  the  dynamics  of  the  media  and  source-receiver 
motions.  This  is  especially  true  when  measurements  are  performed  with  a  moving 
array. 

WAVENUMBER  MEASUREMENTS 

The  signal  coherence  can  be  inferred  from  the  measurement  of  coherent  gain 
realized  by  the  summation  of  n  hydrophones  at  a  given  spatial  wavenumber  k 

N 

S  (w,k)  =  £  x  (w)*A  exp(i<j>  (k))  (61) 

P  n=l  n  n 

where  A^  is  the  array  amplitude  shading  coefficient  and  <f>^(k)  is  a  phase  cor¬ 
rection  to  steer  the  array  in  the  k  direction.  An  example  of  the  beam  output 
|S  (to,k)|^  is  shown  in  figure  14.  In  figure  14a,  one  clearly  sees  the  separa¬ 
tion  of  signal  energy  and  tow-ship  noise  and  in  figure  14b  the  beam  response  of 
the  high  signal  level  173  Hz  tone.  The  estimation  of  received  coherent  signal 
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Figure  13.  Magnitude  squared  coherence  (a)  signal-to-noise  ratio 
(b)  versus  length  along  the  streamer. 
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energy  is  easily  performed  at  these  beam  signal-to-noise  ratios.  In  this  exam¬ 
ple,  the  averaging  interval  used  was  too  long  when  one  considers  the  source- 
receiver  motion.  Nevertheless,  the  peak  signal-levels  provided  a  good  measure 
of  gain  if  the  mean  hydrophone  level  is  known  and  a  functional  form  of  the 
coherence  versus  transverse  distance  can  be  hypothesized.  That  is  to  say,  the 
difference  between  20  Log  (N)  and  the  measured  gain  referred  to  as  the  relative 
array  signal  gain  (RASG)  or  the  differential  array  signal  gain  (DASG)  is  deter¬ 
mined  by  the  degree  of  signal  phase  variability  and  consequently  the  degree  of 
coherence.  The  measurement  uncertainties  are  reduced  due  to  the  large  beam 
signal-to-noise  ratios  and  the  increased  number  of  degrees  of  freedom. 

However,  note  in  figure  14b  the  effect  of  beam  broadening  due  to  the  mul¬ 
tipath  spread  of  off  broad  side  arrivals  is  observed.  Whereas  this  wavenumber 
(angle  domain)  measurement  of  spatial  coherence  has  the  advantage  of  increasing 
signal-to-noise  ratio  it  is  also  susceptible  to  multipath  and  array  deformation 
effects . 


Figure  14.  FFT  beam  response  versus  wavenumber  (FFT  beam  number)  for  the 
signals  from  the  175  Hz  moored  source  (a)  and  the  173  Hz  towed 
source . 


41 


TR  8077 


42  43 

Angle  domain  measurements  made  by  W.B.  Moseley  (NRL)  *  with  a  line 
array  of  hydrophones  positioned  well  within  the  sound  channel  in  deep  wa«.er  are 
shown  in  figures  15  and  16.  Figure  15a  shows  the  average  measured  signal  angu¬ 
lar  pattern  (denoted  by  the  solid  curve)  and  the  reference  response  (denoted  by 
the  thin  curve)  which  would  result  from  a  completely  coherent  plan  wave  imping¬ 
ing  upon  a  completely  straight  array.  As  indicated  in  the  figure,  there  are 
several  major  differences  between  the  measured  and  reference  patterns.  The 
array  signal  gain  (indicated  by  the  measured  value  at  zero  degrees)  is  degraded 
by  about  six  dB  relative  to  the  ideal  reference  value  while  the  three  dB  angu¬ 
lar  width  of  the  measured  signal  is  about  twice  as  wide  as  that  of  the  refer¬ 
ence  pattern.  The  measured  distribution  exhibits  broad  secondary  maxima  at 
angles  different  from  the  source  direction  (these  will  be  shown  to  be  due  to 
multipath  interference).  Generally,  the  signal  energy  is  distributed  across 
angles  instead  of  being  restricted  to  a  narrow  main  lobe  with  low  side  lobe 
structure  as  in  the  reference  case.  The  causes  of  this  increased  spreading  in 
angle  are  threefold--acoust ic  multipath  interference,  array  deformations,  and 
acoustic  scattering  due  to  random  sound  speed  variations  in  the  ocean.  The 
effects  of  each  of  these  mechanisms  were  examined  individually  and  the  results 
compared  with  the  measured  signal  angular  pattern  to  determine  their  relative 
influence . 

Horizontal  array  deformations  were  independently  measured  during  the 
experiment,  and  the  lateral  hydrophone  displacements  were  used  in  computing 
deformed  array  beam  patterns.  Figure  15b  gives  the  comparison  between  the 
average  measured  signal  angular  pattern  and  the  computed  average  beam  pattern 
for  the  deformed  array  (the  thin  curve).  It  is  evident  in  the  figure  that  the 
deformations  of  the  array  have  contributed  significantly  to  the  loss  in  array 
signal  gain  and  to  both  the  observed  spread  of  energy  in  angle  and  the  increase 
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Figure  15,  The  average-measured  beam  response  (dB)  versus  sin0  compared  to 
computed  beam  response  for  (a)  the  case  of  an  ideal  array, 

(b)  the  case  of  a  deformed  array,  (c)  multipath  interference,  and 
(d)  azimuthal  scattering. 


Figure  16.  Comparison  of  measured  and  calculated  beam  angular  response  versus 
sin0  between  ±\/2 d 
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in  the  three  dB  width.  However,  this  mechanism  does  not  account  fully  for  any 
of  these  aspects  of  the  measured  signal  distribution  nor  does  it  produce  the 
broad  secondary  maximum  ("shoulders")  of  the  measured  pattern. 

44  45  46 

The  corrected  parabolic  equation  propagation  model  *  5  using  as  input 

the  measured  sound  speed  profiles  between  source  and  receiver  was  used  to  cal¬ 
culate  the  complex  pressure  field  at  the  location  of  the  receiver.  Then  stand¬ 
ard  beamforming  was  employed  for  a  straight,  predominantly  horizontal  but 
slightly  tilted  array  (tilt  was  determined  from  independent  hydrophone  depth 
measurements)  to  determine  its  angular  response  in  the  presence  of  multipath 
interference.  Figure  15c  shows  the  angular  array  output  due  to  multipath 
interference  (the  thin  curve)  in  comparison  with  the  average  measured  signal 
angular  pattern.  The  comparison  indicates  that  a  substantial  portion  of  the 
signal  gain  degradation  is  due  to  multipath  interference,  that  this  receiver 
has  resolved  part  of  the  upcoming  and  downcoming  energy,  and  that  this  energy 
would  appear  at  the  angles  required  to  explain  the  broad  secondary  maxima  or 
"shoulders"  on  the  measured  signal  angular  response. 

The  prediction  of  azimuthal  scattering  effects  results  from  applying  a 

space-angle  Fourier  transform  to  a  spatial  coherence  function  generated  with 

29-30 

the  Beran  and  McCoy  treatment  of  propagation  in  a  highly  anisotropic  ran¬ 

dom  medium  like  the  ocean,  with  the  environmental  factors  computed  from  inde¬ 
pendent  oceanic  measurements.  Figure  15d  shows  the  theoretical  azimuthal  scat¬ 
tering  effects  (the  thin  curve)  in  comparison  with  the  average  measured  signal 
angular  pattern.  As  is  shown  in  the  figure,  random  medium  scattering  contri¬ 
butes  only  a  small  part  to  the  degradation  in  array  signal  gain  and  the  spread¬ 
ing  of  energy  in  angle.  Here,  its  primary  function  is  an  angular  smoothing 
influence  on  the  effects  due  to  array  deformation  and  multipath  interference. 
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As  was  suggested  earlier,  one  must  look  to  a  combination  of  these  mechan¬ 
isms  to  predict  the  measured  angular  response.  The  predicted  pattern  shown  in 
figure  16  was  synthesized  in  the  following  fashion:  with  the  incorporation  of 
appropriate  phase  delays  to  simulate  the  measured  horizontal  array  deforma¬ 
tions,  beamformer  outputs  were  computed  for  the  predominantly  horizontal  but 
tilted  array  located  in  the  complex  pressure  field  given  for  the  presence  of 
multipath  interference  by  the  corrected  parabolic  equation  model;  then  this 
angular  output  was  convolved  with  the  theoretical  azimuthal  scattering  pattern 
to  arrive  at  the  predicted  pattern.  Figure  16  gives  the  comparison  between  the 
predicted  angular  pattern  (thin  curve)  and  the  average  measured  signal  angular 
pattern  (solid  curve).  There  is  now  an  agreement  between  the  predicted  and 
measured  patterns  to  within  2.5  dB  throughout  almost  all  of  the  angular  domain. 
Thus,  the  combination  of  multipath  interference,  array  deformations,  and  acous¬ 
tic  scattering  accounts  for  the  array  signal  gain  degradation,  the  change  in 
three  dB  angular  width,  the  "shoulders"  in  the  measured  pattern,  and  the  gene¬ 
ral  redistribution  of  energy  in  the  angle  domain. 

47 

TEMPORAL  COHERENCE- -FREQUENCY  DOMAIN  RESULTS 

Thus  far  we  have  discussed  primarily  the  spatial  coherence  problem.  Con¬ 
sider  next,  partial  temporal  coherence  and  spreading  in  the  frequency  domain. 
There  have  been  a  number  of  mechanisms  investigated  as  causes  of  temporal 
coherence  degradation  including:  small  scale  turbulence  and  mixing,  surface 
scattering,  mesoscale  features  like  eddys,  internal  waves,  and  relative  source- 
receiver  motion.  Small  scale  turbulence  and  mixing  may  be  very  important  for 
the  high  frequency,  short  range  problem;  however,  below  1000  Hz  and  for  long 
range  situations,  these  effects  do  not  appear  to  strongly  influence  the  re¬ 
sults.  Surface  scattering  produces  acoustic  sidebands  shifted  about  80-100 
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millihertz  from  the  center  frequency.  On  long  time  scales  of  days  to  months, 
mesoscaie  features  (like  eddys  or  rings)  through  their  generation,  motion,  or 
decay  can  cause  temporal  changes  in  the  received  acoustic  field.  Yet  these 
features  by  themselves  are  not  critical  to  the  coherent  temporal  signal  proces¬ 
sor  design  which  is  most  concerned  with  variations  of  time  scales  from  tens  of 
seconds  to  tens  of  minutes.  For  low  frequency,  long  range,  and  short  time 
scales,  the  temporal  coherence  of  the  received  field  is  dominated  by  two  mecha¬ 
nisms.  One  is  the  internal  wave  (gravity  waves  within  the  ocean)  modulation  of 
the  acoustic  path  structure.  The  other  is  a  variation  of  the  deterministic 
multipath  interference  pattern  at  the  receiver  caused  by  relative  source- 
receiver  motion. 

Two  examples  of  temporal  variations  causing  spreading  in  the  frequency 
domain  are  for  the  fixed  source-fixed  receiver  case  where  ocean  medium  interac¬ 
tions  with  the  propagating  acoustic  field  dominate  the  results.  Measurements 
for  spreading  in  the  frequency  domain  of  an  initially  CW  signal  at  406  Hz 

48 

transmitted  from  Eleuthera  to  Bermuda  were  obtained  by  Clark  and  Kronegold. 

The, spectral  analysis  of  the  phase  and  amplitude  temporal  fluctuations  of  the 
received  signal  showed  substantial  energy  had  been  spread  out  to  frequencies 
about  3  millihertz  shifted  from  the  carrier.  A  processor  with  a  bandwidth  of 
order  3  millihertz  would  accept  the  signal  energy  that  has  been  spread  in  fre¬ 
quency  and  yet  filter  out  extraneous  noise  energy.  Another  Eleuthera  to  Ber- 

49 

muda  experiment  was  reported  by  Stanford.  His  spectra  were  similar  to  those 

-2 

of  the  MIMI  data  in  that  the  slope  of  the  phase  spectra  are  comparable  (f  ) 
and  in  that  substantial  energy  has  been  spread  out  to  about  3  millihertz.  For 
the  Eleuthera  to  Bermuda  propagation  path,  the  highest  occurring  Brunt-Vaisala 
frequency  is  of  order  3  millihertz.  Since  internal  wave  energy  is  not  present 
above  the  highest  Brunt-Vaisala  frequency  and  since  these  spectra  show  that 
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most  of  the  significant  energy  is  contained  within  the  3  millihertz  spread, 
this  is  strongly  suggestive  that  the  internal  wave  mechanism  was  the  dominant 
cause  of  spectral  spreading  for  the  fixed  source-fixed  receiver  case. 

Porter  and  Spindel"^  compared  calculations  with  measurements  to  demon¬ 
strate  when  source-receiver  motions  were  dominant.  Three  cases  were  consid¬ 
ered:  a  fixed  source-receiver,  a  slowly  drifting  receiver,  and  a  rapidly 

drifting  receiver.  The  first  case  was  obtained  for  a  fixed  source-receiver 
geometry.  A  weak  scattering  ray  theory  for  internal  wave  interaction  by  Porter 
and  Spindel  predicts  the  phase  rate  spectra  to  be  a  constant  for  frequencies 
less  than  the  Brunt-Vaisala  frequency  with  a  rapid  fall-off  above  this  fre¬ 
quency.  The  predictions  of  their  theory  were  shown  to  be  in  reasonable  agree¬ 
ment  with  the  measurements.  The  second  case  was  the  results  for  a  slowly 
drifting  sonobuoy  (a  drift  rate  of  0.05  knots)  and  the  third  case  was  for  a 
rapidly  drifting  ship  suspended  hydrophone  (a  drift  rate  of  0.27  knots).  The 
data  were  compared  with  the  phase  rate  spectra  predictions  made  by  Porter  and 
Spindel  using  a  theory  developed  by  Dyson,  Munk,  and  Zetler^^  which  assumes 
that  the  spectral  behavior  is  dominated  by  deep  amplitude  fades  caused  by 
changing  random  multipath  interference.  The  change  from  poor  to  good  agreement 
between  the  Dyson,  Munk,  and  Zetler  model  and  the  data  as  the  drift  rate  in¬ 
creases  was  attributed  to  the  increasing  dominance  of  the  source-receiver 
motion  induced  multipath  interference  effects. 

The  precise  magnitude  of  source-receiver  motion  required  to  dominate  a 

particular  temporal  spectrum  has  yet  to  be  completely  specified.  Dashen, 

52 

Flatte,  Munk,  and  Zachariasen  have  estimated  that  radial  source  motion  of 

0.33  knots  and  lateral  motion  of  2.2  knots  will  significantly  influence  the 
53 

spectra.  Dyer  has  estimated  that  source-receiver  motions  of  3  to  5  knots 
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will  produce  dominance  by  multipath  interference  effects.  Results  of  a  treat- 
mend  of  source-receiver  motion  induceu  deterministic  multipath  effects  by 
Clark,  Flanagan,  and  Weinberg  shows  a  spread  in  the  frequency  domain  of  about 
3  millihertz  occurs  for  speeds  of  5  knots.  This  implies  that  at  5  knots  varia¬ 
tions  due  to  source  motion  may  be  comparable  with  those  caused  by  internal  wave 
effects  (assuming)  that  the  source  and  receiver  depths  and  sound  speed  profile 
are  appropriate).  As  the  speed  increases,  the  frequency  spreading  increases  so 
that  for  high  rates  of  source  motion  a  very  narrow  bandwidth  processor  would 
filter  out  some  of  the  source  energy.  Thus,  both  internal  waves  and  source- 
receiver  motion  can  limit  the  useful  frequency  resolution  of  the  processor. 

In  summary,  the  topics  of  acoustic  spatial  coherence  and  acoustic  temporal 
coherence  have  been  treated  in  a  conceptually  simple  but  analogous  fashion  to 
indicate  that  they  are  concerned  with  the  spreading  of  the  acoustic  energy  in 
angle  and  frequency,  respectively.  For  the  long  range,  low  frequency,  short 
time  scale,  predominantly  refracted  and  refracted-surface-ref lected  energy 
reception  situation,  the  physical  processes  and  environmental  influences  which 
can  degrade  coherence  have  been  indicated  when  the  spatial  and  temporal  effects 
are  separable.  Three  factors  for  the  spatial  case--acoustic  multipath  inter¬ 
ference,  array  deformations,  and  scattering  in  a  random  medium--and  two  factors 
for  the  single  hydrophone  temporal  case--internal  waves  and  source-receiver 
relative  motion--were  shown  to  directly  limit  the  resolution  capability  in 
angle  and  frequency  and  thereby  the  receiving  system  gains.  Over  the  past  few 
years  theories  and  experiments  have  contributed  to  a  better  understanding  of 
the  phenomena  and  have  given  verification  in  uncomplicated  situations.  How¬ 
ever,  usually  the  effects  of  the  various  factors  are  modeled  separately.  Theo¬ 
ries  need  to  be  developed  both  for  spatial  and  temporal  coherence  which  treat 
several  factors  simultaneously. 
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THE  OBSERVATION  OF  MULTIPATHS 


The  propagation  of  mid-to-low  frequency  sound  in  the  ocean  to  long  range 
with  refracted-refracted  (RR)  or  refracted-surface-ref lected  (RSR)  paths  has 
been  observed  by  several  investigators  (55,56)  with  arrays  which  separated  the 
multipath  arrivals.  The  separability  of  these  arrivals  by  either  a  vertical  or 
horizontal  array  is  determined  by  the  length  of  the  array,  its  orientation,  and 
the  coherence  length  of  the  signal.  The  effective  coherence  length  of  the  sig¬ 
nal  is  determined,  in  this  instance,  by  the  amount  of  angular  spread  in  the 
signal  caused  by  diffused  energy  about  a  single  arrival  direction  or  by  the 
distribution  of  energy  in  angle  caused  by  the  addition  of  several  paths.  Con¬ 
sider  a  plane  wave  from  the  0  ,<j>  direction  with  wavenumber  vector  k  incident 

o  o  — 

on  an  array  along  the  x  axis  with  a  slight  downward  tilt  as  shown  in  figure 
17. 

k*x/|k||x|  =  cosG  =  cos0  cos<i>  cosiIj  +  sin0  sinii;  =  sin0’  (62) 

- 1  1  1  a  o  o  o  a 


cos  da  =  cos (p0  cos 60  cos ijj  +  sin#0  sin<// 

Figure  17.  The  geometry. 
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We  have^  shown  from  the  of f -broadside  response  of  a  uniform  line  array 

2AsinG '  =  2Acos6  =  2A/L  (63) 

a  a 

between  the  first  nulls  and  for  the  half  power  (hp)  points 

(AsinG')  =  (AcosG)  *  0.886A/L  (64) 

a  hp  hp 

1/2 

and  at  end-fire  (AG).  =  2(0.886A/L)  .  For  the  off-broadside  case,  one  can 

hp 

calculate  a  matched  field  condition  where  the  angular  spread  of  the  signal  is 
equal  to  the  3  dB  beamwidth  of  the  line  array  when  the  beam  is  steered  in  the 
direction  of  the  signal. 

L/(A/AcosG  )  =  Lcos(f>  cos^/(A/AcosG  )  +  Lsin^/( A/AsinG  )  (65) 

a  o  o  o 

The  quantities  A/A  cos  0  88  L  A/AcosG  =  L  ,  A/AsinG  =  L  are  inter- 

o  eff  ox  o  z 

preted  as  coherence  lengths.  When  <  LcosG^  cost/;,  the  horizontal  component 
of  the  array  over-resolves  the  field  as  the  beam  width  for  a  given  conical  an¬ 
gle  is  narrower  than  the  field.  When  LsiniJ;  >  A/Asin0  ,  the  vertical  resolution 

o 

of  the  array  due  to  its  vertical  extent  over-resolves  the  field.  When  equation 
65  is  satisfied  one  has  matched  the  field  to  the  length  of  the  system. 

The  angular  spreads  could  be  due  to  the  distribution  of  energy  about  each 
path  resulting  from  scattering,  the  different  vertical  arrival  angles  of  inde¬ 
pendent  paths,  or  the  angular  spread  across  the  array  due  to  a  change  in  range. 

The  case  of  vertical  arrival  of  different  paths  in  the  deep  ocean  can  be 

57  58 

estimated  by  use  of  Snell's  law  or  Cox's  ARAD  5  technique.  Calculations 

with  a  representative  Atlantic  profile  (V  «  1537  m/sec,  V  =  1526  m/sec,  V,  = 

os  a  b 

1554  m/sec)  and  Snell's  Law  yield  for  the  bottom  grazing  and  surface  grazing 
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rays  Lx  =  90,9  X  and  =  14,4  X.  Using  the  ARAD  technique  with  AcosO  =  V/V  , 

A/D2  =  0.17,  and  D(km)  depth  excess  results  in  L  =89  A/D  and  L  =18  A/d. 

x  z 

The  deterministic  multipath  effects  can  be  described  by  the  methods  of 
geometrical  acoustics  at  the  higher  frequencies  and  by  modes  at  the  lower  fre¬ 
quencies,  The  equivalence  of  modes  and  rays  has  been  discussed  by  Guthrie  and 
Tindle.“^  Here  we  use  the  ray- theoretic  model,  Trimain^  and  the  normal  model 
mode  of  Bucker^  to  illustrate  these  effects. 

Figure  18a  shows  the  path  arrival  structure  as  computed  with  the  Trimain 
Code^  at  a  frequency  of  175  Hz  for  a  source  depth  of  430  m  and  a  receiver 
depth  of  250  m.  The  ordinate  is  range  in  0.2  km  increments  and  the  abcissa  is 
the  vertical  arrival  angle  with  0°  being  the  horizontal.  At  each  range  incre¬ 
ment,  the  vertical  line  length  represents  the  intensity  for  the  eigenray  at  the 
angle  corresponding  to  the  intersection  of  the  vertical  line  with  the  horizon¬ 
tal  range  line.  Observe  the  angular  difference  and  the  convergence  zone 
effects.  The  rate  of  change  of  arrival  angle  with  range  shows  a  definitive 
pattern  with  slopes  between  0.44°/km  to  0.16°/km,  Shown  in  figure  18b  is  a 
similar  case  for  the  Ionian  Basin  of  the  Mediterranean.  In  this  case  the 
source  is  at  a  depth  of  300  m,  the  receiver  is  at  a  depth  of  250  m  and  the  fre¬ 
quency  is  175  Hz.  This  pattern  is  also  distinctive  with  clear  separable  arri¬ 
vals.  The  response  of  a  slightly  tilted  horizontal  array  to  the  field  is  the 
convolution  of  the  beam  response  pattern  and  the  angular  spectrum  of  the  inci¬ 
dent  acoustic  field.  At  any  given  range  one  can  observe  multiple  arrivals  the 
structure  of  which  changes  as  the  source-receiver  separation  is  increased  or 
decreased . 
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Fig.  18a.  Computed  eigen  ray  arrival 
angle  and  intensity  as  a 
function  of  range,  Atlantic 
Basin  at  175  Hz. 


Fig.  18b.  Computed  eigen  ray  arrival 
angle  and  intensity  as  a 
function  of  range,  Ionian 
Basin  at  175  Hz. 


Figure  19  shows  measured  data  obtained  by  Williams  and  Fisher  compared 
to  simulated  data.  The  data  are  for  the  Pacific  with  a  near  surface  source  and 
a  receiver  in  the  upper  part  of  the  sound  channel  in  the  mid-frequency  range. 
The  arrival  angles  range  between  8°  and  14°  in  a  clearly  distinctive  fashion 
with  a  range  rate  of  0.143°/km.  Thus  the  picture  which  data  and  calculation 
show  is  one  of  multipaths  at  different  angles  which  have  definitive  range  rate 
patterns.  Similar  results  have  been  obtained  by  Lawrence  and  Ramsdale.^^ 

These  multipath  patterns  have  been  observed  to  change  with  time  as  the 
properties  of  the  intervening  water  column  changes.  These  temporal  effects 
have  been  discussed  in  the  preceeding  sections.  The  question  of  the  response 
of  an  array  moving  in  these  fields  is  discussed  later. 
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Figure  19. 


Measured  and  computed  arrival  angle  versus  range.  (a)  Measured 
arrival  angles.  (b)  Calculated  arrival  angles  versus  range. 
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SYNTHETIC  APERTURE  RESULT62 

Experimental  results,  shown  in  figure  20,  demonstrate  the  formation  of 
synthetic  apertures  when  the  source-receiver  motion  is  known  or  inferred  by  a 
velocity  filter,  and  when  the  synthetic  aperture  length  is  less  than  the  signal 
coherence  length,  and  when  the  processing  time  is  less  than  the  signal  decorre¬ 
lation  time. 

Synthetic  apertures  have  been  used  in  the  field  of  radar  since  1953.  A 

63 

monograph  of  classic  papers  in  radar  has  been  published  by  Kovaly.  Early 

attempts  at  the  formation  of  synthetics  in  underwater  acoustics  concentrated  on 

kHz  frequencies  and  active  systems.  Sam  Hanish  (NRL)  reviewed  this  field  in 

64 

1975  and  Louis  J.  Cutrona  published  a  paper  in  August  of  1975  which  discussed 

39 

active  sonar  system  design.  In  1976,  observations  by  Ross  Williams  and  R. 

40 

Fitzgerald  .  based  on  single  element  phase  records  led  to  the  speculation  and 
conclusions  that  passive  synthetic  apertures  could  be  formed.  The  Fitzgerald 
letter  focused  on  a  10  Hz  (150  m  wavelength)  signal  over  a  distance  of  65  km 
(R/A  =  6.5  x  103)  and  the  inference  that  apertures  up  to  3300  m  (L/A  =  22) 
could  be  formed.  This  conclusion  was  not  demonstrated  by  coherent  summation  of 
the  signals  or  by  the  measurement  of  signal  gain. 

Ross  Williams  discussed  phase  versus  time  records  obtained  by  a  moored 
sensor  and  a  moving  source  for  distances  up  to  508  km  (R/A  =  104-105)  at  a  fre¬ 
quency  of  400  Hz.  He  was  able  to  find  a  linear  trend  between  the  average  phase 
and  time  for  several  intervals  during  which  he  speculated  that  apertures  up  to 
900  m  (240  L/A)  could  be  formed.  The  effects  of  phase  fluctuation  and  sensi¬ 
tivity  of  the  technique  to  source-receiver  motion  were  discussed  by  the  author 
but  were  not  incorporated  into  an  estimate  of  coherent  gain. 
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(a) 


(b) 


AZIMUTHAL  ANGLE 


Figure  20.  Results  of  a  synthetic  aperture  formed  from  sub-apertures  with 
four  hydrophone  groups. 


55 


TR  8077 


Spindel  and  Porter  have  also  reported  measurements  at  two  frequencies 
(220-406  Hz)  with  drifting  hydrophones  over  a  source-receiver  range  of  260  km 
and  showed  during  periods  of  highly  correlated  phase  the  formation  of  synthetic 
aperture  length  up  to  8  km. 

The  data  shown  in  figure  20  were  obtained  with  four  hydrophone  groups 

6  2 

spaced  0.3  X  apart.  The  response  of  the  subaperture  array  was  determined  by 

N 

b(k  , 0  } t  )  =  Y  S  exp(i(cot  -  k*r  ))  A  exp( ikirdncosS  )  (66) 

sso  %or  o - nnr  s 

n=l 

represents  the  shading  coefficient  taken  as  unit  and  0^  is  the  phase 
steering  angle  which  determines  the  azimuthal  response  of  the  array.  The  for¬ 
mation  of  the  synthetic  aperture  is  accomplished  by  adjusting  the  successive 
beamformed  outputs  by  a  factor  $  which  should  equal  jc-VmAt  in  the  direction  of 
the  arrival.  The  synthetic  aperture  result  can  be  shown  to  be 

M=1 

B(k  ,t  ,t  +nAT)  =  [b(k  ,0  ,t  )]  Y  exp(i(a)mAt  -  k*VmAt  -<(>))  (67) 

soo  sso  u  ^  -  m 

m=0 

where  mAt  is  the  number  of  time  intervals,  k_  is  the  wavenumber  vector  and  v_  is 
the  source  receiver  velocity.  The  term  unnAt  accounts  for  the  phase  due  to  the 
time  difference  between  t^  and  t  ^ .  The  term  _k*VmAt  accounts  for  the  spatial 
distance  VmAt  =  AL  and  the  direction  wavenumber  vector  k. 

The  results  shown  in  figure  20a  show  the  response  of  the  1.2  X  array,  for 
which  the  arrival  structure  is  not  resolved.  Second,  figure  20b  is  the  23  X 
array  which  shows  the  direction  of  the  source  and  an  aliased  side  lobe  due  to 
the  2  X  subaperture  spacing.  Figure  20c  shows  the  95  X  array  which  begins  to 
resolve  the  arrival  structure.  The  aliased  response  is  also  evident.  These 
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data  show  when  the  single  path  coherence  length  is  greater  than  the  length  of 
the  array  and  the  temporal  coherence  length  is  greater  than  the  processing  time 
that  synthetic  apertures  can  be  formed.  The  result  shows  the  separation  of 
arrival  structure  as  the  resolution  of  the  array  is  increased.  The  off-broad¬ 
side  arrivals  are  seen  to  be  affected  by  the  vertical  arrival  structure.  This 
result  strongly  implies  that  motion  between  the  source  receiver  can  also  affect 
the  beam  response  output. 

DYNAMIC  EFFECTS 

The  space-time  measurement  of  the  ocean  acoustic  pressure  field  has  been 
described  primarily  for  stationary  sources  and  receivers.  In  these  cases,  the 
limiting  resolution  of  the  line  array  was  shown  to  be  due  to  the  angular  spread 
of  the  incident  energy  about  a  single  arrival  due  to  scattering  phenomena  or 
due  to  the  addition  of  multipaths  and  their  variation  over  the  aperture.  Pre¬ 
viously,  evidence  obtained  by  Spindel  and  Porter  as  well  as  estimates  by  Dyer 
on  the  importance  of  multipath  fluctuations  and  its  dominance  as  relative 
motion  between  the  source  and  receiver  exceeds  1.5  m/sec  was  described.  In 
these  instances,  the  temporal  and  spatial  resolution  relationships  are 


Af-AT  ^  N 

(68) 

AG  sinG  (L/A)  =  AG*  cosG'CL/A)  ^  1 
a  a  a  a 

(69) 

where 

sinG  =  cosG '  and  cosG  =  cosd>  cosG  (70) 

a  a  a  o  o 

These  relationships  follow  directly  from  our  use  of  the  Fourier  transform 
and  a  rectangular  sampling  window.  We  observe  that  these  relationships  can 
become  coupled  by  motion  during  the  observation  time. 
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A0 

oa 


90 


oa 

9R 


AR  + 


90 


oa 

9<j> 


A<J> 


sin0  cos  d> 
_ o _  o 

sin0 

a 


90 

if  <L<:os+0  +  VT) 


sin0  cos0 

+  Sine  °  (VT/E)  (71) 

a 

This  relationship  can  also  be  applied  to  the  multipath  problem  as  long  as 

the  angular  separation  between  the  paths  is  not  great.  We  let  A0  represent 

oa 

the  angular  spread  due  to  the  environmental  factors  and  A0^  the  beamwidth  of 

the  measurement  system.  When  these  two  quantities  are  equal,  we  have  matched 

the  measurement  system  to  the  environment.  When  A0  <  A0  .  we  have  the  capa- 

a  oa 

bility  to  resolve  the  multipaths  and  when  the  acoustic  intensity  is  contained 
within  a  beam.  A  measurement  problem  of  considerable  interest  is  the  resolu¬ 
tion  of  distant  shipping.  In  the  following  we  will  assume  radial  (V  )  and  azi- 

R 

muthal  velocities  (V^)  of  10  m/sec,  an  azimuthal  bearing  of  =  80°,  90°  being 
broadside,  and  measurement  interval  of  300  seconds. 

By  equating  the  aperture  width  with  the  angular  spread,  it  can  be  shown 
that  for  a  single  path  with  a  radial  motion  f 

X/L  «  C  g  (0  ,  <f>  )  (Lcos<f>  +  V  AT)  (72) 

e°o  o  o  o  R 

8o(0o>  ♦«,>  =  sin6o  cosV  ce  =  l"o/9Rl 

The  constant  is  the  rate  of  change  of  single  path  vertical  arrival 
angle  (0Q)  and  can  be  taken  from  the  previously  shown  vertical  arrival  angle 
plots.  This  constant  is  small  0(2x10  ^)  and  near  convergent  zones,  0(8x10  ~>). 

The  azimuthal  variation  for  the  single  path  case  can  also  be  shown  to  be: 

X/L  =  sin<f)  cos0  V  AT/R  (73) 

o  o  <p 


58 


TR  8077 


Where  R  is  the  radial  distance  to  the  source  taken  here  to  be  70  km.  Sim¬ 
ilar  results  can  be  obtained  for  the  multipath  case  and  are  given  by 

3A0  30 

X/L  =  cosO,  AcosO  {1  +  (— )  |-^|  +  VrAT} 

O 

cos0  sin<() 

+  - ^ - -  V  AT  (74) 

R  < p 

g  =  tan2(J>  /  tan0  n  (tan2<f>  +  sin20  ,) 

1  o  ol  o  ol 

The  multipath  case  shows  a  constant  term  determined  by  the  initial  spread 
in  angle  between  the  paths,  a  radial  term  which  accounts  for  the  change  in 
spread  with  range,  a  radial  term  which  accounts  for  the  variation  in  central 
arrival  angle  with  range  and  finally  a  term  which  accounts  for  the  azimuthal 
rate  of  change. 

These  results  are  interesting  as  they  allow  us  to  rank  the  relative  impor¬ 
tance  of  these  effects  on  the  ocean  measurement  problem.  When  motion  is  not 
important  the  single  path  refraction  result  yields 

L/X  +  8  x  103//T 

Comparison  with  results  from  the  “Experimental  Results"  section  indicates  that 
in  this  instance  the  measurement  will  be  dominated  by  volume  or  surface  scat¬ 
tering.  When  the  motions  are  extreme  we  find 

L/X  (1/C  )  g  VAT  =  3-104 
e  o  R 

and  again  scattering  not  refraction  limits  the  ability  to  measure  the  single 
path  arrival.  These  results  confirm  our  previous  discussion  regarding  single 
path  propagation  and  volume  scattering. 
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We  can  continue  the  process  and  estimate  term  by  term  the  relative  impor¬ 
tance  of  each  effect.  For  the  specific  example  given  here,  we  have  tabulated 
these  results  in  the  following  table. 


A  Comparison  of  Multipath/Motion  Effects 


CASE 

L/A 

REMARKS 

Single  Path 

Static 

8x103//a" 

*  0  Scattering  is 

dominant . 

Dynamic  (radial) 

4xl03 

V  >  20 
r 

Dynamic  (azimuthal) 

25 

at  R  =  70  km,  azimuthal 
very  important 

motion  is 

Multipath 

Static  (Atlantic) 

525 

o 

O 

00 

II 

o 

-0- 

105 

<f>  =  30° 

o 

506/D  =  84 

Cox's  ARAD,  D(km) 

Dynamic  V  =  V,  =  0 
r  <f> 

V  =  0,  V  =20 
<p  r 

302 

278 

Here  we  use  a  case  with 
602  "  16° 

0O1  "  8'5°’ 

V  =  0,  V  =20 
<p  r 

67.3 

CZ  case 

o 

ii 

u 

> 

o 

CN 

II 

> 

.  21.6 

The  results,  although  only  given  for  a  specific  measurement  problem, 
clearly  show  the  importance  multipath  arrival  structure  and  motion  have  on  the 
resolution  capabilities  of  the  ocean  measurement  system.  Measurements  per¬ 
formed  with  arrays  to  measure  horizontal  noise  directionality  should  carefully 
consider  these  results.  We  have  used  only  horizontal  arrays  with  no  tilt  or 
vertical  extent.  These  results  can  easily  be  extrapolated  to  the  case  of  a 
tilted  array  and  the  procedure  applied  to  the  case  of  a  vertical  array.  The 
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results  thus  far  have  used  the  vertical  arrival  structure  for  the  Atlantic, 
figure  18a.  These  results  scale  to  the  cases  of  the  Pacific,  figure  19,  and 
the  Mediterranean,  figure  18b. 

Figure  21  illustrates  a  numerical  calculation  using  a  normal  mode  model 
and  a  Mediterranean  sea  profile.  The  array  has  a  slight  tilt  and  the  beam 
response  and  illumination  function  are  plotted  versus  range.  We  observe  the 
variation  with  range  in  the  array  illumination  pattern  due  to  the  combinations 
of  multipath  and  the  subsequent  bifurcation  of  the  beam  response.  The  results 
are  plotted  in  1  km  intervals  and  provide  a  graphic  illustration  as  to  varia¬ 
tion  expected  in  range  as  a  function  of  source  receiver  motion.  The  range 
change  during  a  measurement  period  of  300  seconds  is  154.2  m/knot.  Since  the 
variations  with  range  are  on  the  order  of  1  km,  we  can  expect  dramatic  varia¬ 
tions  as  our  speed  approaches  3  knots.  This  result  is  wholly  consistent  with 
our  previous  comments  regarding  multipath  fluctuations. 

The  conclusion  drawn  here  is  that  the  resolution  of  multipath  arrivals 
with  a  horizontal  array  with  slight  tilt  occurs  well  before  single  path  limits 
are  reached.  Furthermore,  relative  motions  of  the  source  of  sound,  i.e.,  a 
surface  ship,  further  impact  the  choice  of  measurement  parameters  available  at 
a  given  frequency,  that  is,  the  array  length,  and  the  observation  time.  These 
results  should  apply  to  many  ocean  measurement  problems  such  as  the  determina¬ 
tion  of  the  role  of  shipping  on  the  directionality  of  noise. 
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ILLUMINATION  OF  APERTURE 


Figure  21.  Calculated  near-broadside  beam  response  patterns  and  illumination 
of  aperture  as  a  function  of  range. 
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SUMMARY  AND  CONCLUSIONS 

This  paper  had  as  its  focus  the  problem  of  ocean-acoustic  measurements 
with  line  arrays.  The  space- time  transforms  were  reviewed  to  demonstrate  the 
analogy  between  spatial  and  temporal  properties,  to  stress  the  importance  of 
convolution  and  matched  field  processing.  A  criteria  was  presented  by  which 
the  resolution  of  such  measurement  systems  could  be  calculated. 

The  static  source-receiver  case  was  shown  to  be  influenced  by  the  random¬ 
ness  in  signal  phase  due  to  scattering  and  closely  following  the  method  and 
approach  of  Shifrin,  we  show  the  results  of  this  Gaussian  phase  randomness  and 
the  importance  of  the  coherence  functional  form.  Calculations  and  data  were 
employed  to  show  the  importance  of  multipath  effects  on  the  relative  gain  of 
line  array  measurements  systems  and  the  difficulties  encountered  for  the  deter¬ 
mination  of  coherence  lengths.  Single  path  coherence  lengths  were  found  to  be 
large  and  predictable  using  an  environmental  parameter  E  and  the  Beran-McCoy 
mutual  coherence  function  of  form.  Nevertheless,  multipath  effects  were  in  our 
opinion  dominant. 

The  temporal  fluctuation  problem  was  briefly  introduced  for  completeness 
but  also  to  stress  the  finding  that  for  relative  source-receiver  speeds  of  1.5 
m/sec  (3  knots)  or  greater,  the  fluctuations  are  dominated  by  the  changes  in 
the  multipath  arrivals.  These  effects  dominate  most  measurements.  We  illus¬ 
trate  slowly  moving  source-receiver  and  show  with  data  that  a  line  array  can 
resolve  distinct  multipath  arrivals.  We  presented  a  simplified  analytical  for¬ 
mulism  for  the  determination  of  the  aperture  length  required  to  contain  the 
multipath  arrival  during  the  measurement  period  or  for  the  determination  of  the 
length  of  the  system  required  to  resolve  two  distinct  arrivals. 
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In  closing,  we  conclude  that  vertical  and  horizontal  line  arrays  are  mere¬ 
ly  spatial  transforms  analogous  to  temporal  transforms.  The  spatial  structure, 
angular  spread,  determines  the  length  of  the  antenna  required  either  to  charac¬ 
terize  the  structure  or  to  contain  the  structure  within  a  spatial  wavenumber 
interval . 

When  the  incident  signal  is  a  wave  with  a  normal  variable  in  phase,  the 
net  gain  was  shown  to  proceed  according  to  a  2L/X  until  the  spread  in  incident 
energy  becomes  greater  than  the  spatial  beamwidth  half  power  points,  that  is,  a 
saturation  effect.  We  conclude  that  this  beginning  of  saturation  is  one  limit 
on  measurement  system  resolution. 

When  the  vertical  arrival  angles  from  multiple  paths  project  horizontal 
angles  greater  than  the  half-power  beam  widths,  beam  bifurcation  occurs.  For 
vertical  arrays,  we  simply  have  the  arrivals  on  different  beams.  Horizontal 
arrays  have  an  azimuthal  variation  in  the  resolution  from  broadside  to  end 
fire.  Our  calculations  and  data  result  in  the  conclusion  that  the  combined 
effect  of  antenna  distortion  and  multipath  arrival  structure  are  exceedingly 
important  in  the  ability  of  the  system  to  measure  the  incident  field. 

Motion  of  the  source-receiver  with  respect  to  one  another  couples  the 
space-time  processing  problem.  We  concluded  that  for  static  cases  we  must 
restrict  our  processing  to  temporal  and  spatial  coherence  lengths  determined  by 
either  the  temporal  or  spatial  correlation  function.  However,  when  the  rela¬ 
tive  motion  becomes  important  (>1.5  m/sec),  we  are  restricted  to  an  effective 
length  determined  by  the  time  required  to  perform  the  measurement  and  the  vari¬ 
ation  of  the  multipath  arrival  structure. 
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The  measurement  of  the  signal  gain  of  a  spatial  array  requires  either  the 
determination  of  the  phase  moments,  the  spatial  coherence  or  the  coherent  out¬ 
put  of  the  line  array.  Whereas  it  is  difficult  to  measure  signal  phase  or 
coherence  at  low  signal-to-noise  ratios  and  large  separations,  it  is  possible 
to  measure  the  coherently  summed  output  of  the  array.  This  sum  can  be  related 
to  the  mutual  coherence  function  and  consequently  the  coherence  length  esti¬ 
mated.  Thus,  high  resolution,  large  arrays  may  provide  a  tool  by  which  spatial 
coherence  lengths  can  be  estimated  in  the  low-  to  mid-frequency  region. 

In  the  case  of  the  spatial  coherence  problem,  a  theory  is  required  that 
combines  and  includes  both  the  scattering  and  the  multipath  interference 
effects  simultaneously.  For  the  temporal  coherence  problem,  a  theory  is  needed 
that  incorporates  both  the  internal  wave  modulation  of  the  acoustic  path  struc¬ 
ture  and  the  variation  of  the  multipath  interference  pattern  due  to  relative 
source-receiver  motion.  These  are  current  areas  of  activity  by  several 
research  groups.  A  complete  theory  which  treats  both  temporal  and  spatial 
coherence  simultaneously  (including  the  several  influential  factors  for  each) 
is  in  all  likelihood  some  years  away  from  development.  However,  a  precise  pre¬ 
dictive  capability  for  either  or  both  of  the  acoustic  field  coherence  proper¬ 
ties  will  influence  both  array  and  processor  usage. 
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